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Abstract 

We present a calculation of next-to-leading-order (NLO) QCD corrections to total hadronic 
production cross sections and to light-hadron-decay rates of heavy quarkonium states. Both 
colour-singlet and colour-octet contributions are included. We discuss in detail the use of 
covariant projectors in dimensional regularization, the structure of soft-gluon emission and 
the overall finiteness of radiative corrections. We compare our approach with the NLO ver- 
sion of the threshold-expansion technique recently introduced by Braaten and Chen. Most 
of the results presented here are new. Others represent the first independent reevaluation of 
calculations already known in the literature. In this case a comparison with previous findings 
is reported. 
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1 Introduction 



Since the discovery of the J/ip [Q and its interpretation within QCD as a charm-anticharm bound 
state 1^, the study of quarkonium has received much attention from both a theoretical and an 
experimental point of view, providing a good testing ground for studies of Quantum Chromody- 
namics (QCD) in both its perturbative and non-perturbative regimes. 

Decay rates of heavy quarkonium states into photons and light hadrons have first been calcu- 
lated at leading-order (LO) and compared to experimental data (see, e.g., [§]) under the assump- 
tion of a factorization between a short- distance part describing the annihilation of the heavy-quark 
pair in a colour- singlet state and a non-perturbative long-distance factor, related to the value at 
the origin of the non-relativistic wave-function or its derivatives. Calculations of decay rates at 
full next-to-leading order (NLO) have also been performed in this framework since the early days 
of quarkonium physics [§]. This approach can be extended to the case of quarkonium production. 
The prescription to evaluate the short-distance coefficients is very simple: the QQ pair has to be 
produced during the short-distance interaction in a colour-singlet state, with the same spin and 
angular momentum quantum numbers of the quarkonium state we are interested in 0. A single 
non-perturbative parameter, the same appearing in the quarkonium decay and provided by the 
bound state Bethe-Salpeter wave-function, accounts for the hadronization of the QQ pair into the 
physical quarkonium state. Applications of this approach, usually referred to as the Colour Sin- 
glet Model (CSM), led to the calculation of LO matrix elements for the production of total cross 
sections and pt distributions in hadronic collisions [^]. Phenomenological successes and failures 
of the CSM in describing the data available up to 1993 on charmonium production and decay are 
nicely reviewed in ref. [0]. 

In the case of S'-waves the simple factorization hypothesis underlying the CSM was confirmed 
by the calculation of one loop corrections ^. The appearance of a logarithmic infrared diver- 
gence in the case of NLO P-wave decays into light hadrons |10| violated instead factorization 
explicitly. Phenomenologically, this singularity could be handled by relating it to the binding en- 
ergy of the bound quarks. Nevertheless, it has since then been clear that in spite of its simplicity 
and physical transparency, the CSM suffers from serious theoretical limitations. The principal 
one being the absence of a general theorem assessing its validity in higher orders of perturbation 
theory, as already indicated by the explicit example of P-states decay. The striking observation 
by CDF of large-pT J/ip and ip' states produced at the Tevatron [jTT], |12| at a rate more than one 
order of magnitude larger than the theoretical prediction, and the serious discrepancies between 
fixed-target data and predictions for the relative production rates of ^Si ^Pi and ^P2 states (for 
recent total cross-section measurements, see |]13|, |14|); provided equally compelling evidence that 
some important piece of physics was missing from the CSM. 

The road to solve these formal and phenomenological problems has been indicated by the work 
of Bodwin, Braaten and Lepage (BBL) ||T^, which provided a new framework for the study of 
quarkonium production and decay within QCD. In this work, perturbative factorization is retained 
by allowing the quarkonium production and decay to take place via intermediate QQ states with 
quantum numbers different than those of the physical quarkonium state which is being produced 
or which is decaying. In the case of production, for example, the general expression for a cross 
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section is given by: 

da[H + X)=Y,da{QQ[n]+X){0''[n]) . (1) 

n 

Here da{QQ[n] + X) describes the short distance production of a QQ pair in the colour, spin 
and angular momentum state n, and {0^[n\), the vacuum expectation value of a four-fermion 
operator defined within Non-Relativistic QCD (NRQCD) [^6], describes the hadronization of 
the pair into the observable quarkonium state H. QQ states with quantum numbers other than H 
arise from the expansion of the H Fock-space wave function in powers of the heavy quark velocity 
V. The relative importance of the various contributions in eq. (|l|) can be estimated by using 
NRQCD velocity scaling rules , which allow the truncation of the series in eq. (0) at any given 



order of accuracy. If one only retains the lowest order in v, the description of S'-wave quarkonia 
production or annihilation reduces to the CSM one. In the case of P waves, instead, contributions 
from colour-octet, S'-wave QQ states are of the same order in v as those from the leading colour- 
singlet P-wave state. Infrared singularities which appear in some of the short-distance coefficients 



of P-wave states can then be shown |]T5|] to be absorbed into the long-distance part of colour-octet 
S'-wave terms, thereby ensuring a well defined overall result. This framework for the calculation of 
quarkonium production and decay is often referred to as the colour-octet model (COM), perhaps 
with an abuse of language, as the COM pretends to be a direct outcome of QCD, rather than just 
a "model". Good reviews of the underlying physical principles and applications can be found in 



refs. |T8|, m 



The effect of colour-octet contributions can be extremely important even in the case of S'-wave 
production. In fact, while their effects are predicted by the scaling rules to be suppressed by 
powers of v with respect to the leading colour- singlet ones, their short-distance coefficients can 
receive contributions at lower orders of a^, thereby enhancing significantly the overall production 



rate |20|. The inclusion of these colour-octet processes, in conjunction with the observation that 



gluon fragmentation provides the dominant contribution to the short-distance coefficients at large- 



Pt leads to a very satisfactory description of the Tevatron data [^, 

One of the most important consequences of factorization for quarkonium production is the 
prediction that the value of the non-perturbative parameters does not depend on the details of the 
hard process, so that parameters extracted from a given experiment can be used in different ones. 
For simplicity, we will refer to this concept as "universality" . Several studies of experimental data 
coming from different kind of reactions have been performed to assess the validity of universality. 
For example calculations of inclusive quarkonia production in e^e~ fixed target experiments 
26| , 7p collisions [^, ^ and B decays p9| have been carried out within this framework. The 



overall agreement of theory and data is satisfactory, but there are clear indications that large 
uncertainties are present. The most obvious one is the discrepancy |^ between HERA data 
and the large amount of inelastic J/ip photoproduction predicted by applying the colour-octet 
matrix elements extracted from the Tevatron large-p^ data ||2^, |3T1] . 

It becomes therefore important to assess to which extent is universality applicable. Several 
potential sources of universality violation are indeed present, both at the perturbative and non- 
perturbative level. On one hand there are potentially large corrections to the factorization theorem 
itself. In the case of charmonium production, for example, the mass of the heavy quark is small 
enough that non-universal power-suppressed corrections can be large. Furthermore, some higher- 
order corrections in the velocity expansion are strongly enhanced at the edge of phase-space . 
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For example, the alternative choices of using as a mass parameter for the matrix elements and for 
the phase-space boundary the mass of a given quarkonium state or twice the heavy-quark mass 
2m, give rise to a large uncertainty in the production rate near threshold. These effects, which are 
present both in the total cross-section and in the production at large-pr via gluon fragmentation, 
violate universality. This is because the threshold behaviour depends on the nature of the hard 
process under consideration. 

Another source of bias in the use of universality comes from purely perturbative corrections. 
Most of the current predictions for quarkonium production are based on the use of leading- 
order (LO) matrix elements. Possible perturbative i^-factors are therefore absorbed into the 
non-perturbative matrix elements extracted from the comparison of data with theory. Since the 
size of the perturbative corrections varies from one process to the other, an artificial violation of 
universality is introduced. Examples of the size of these corrections are given by the large impact 
of fcy-kick effects and initial-state multiple-gluon emission in open-charm [^] and charmonium 
production 



In this paper we focus on the evaluation of the 0{al) corrections to quarkonium total hadro- 
production cross sections. Of all the efforts needed to improve the quantitative estimates of 
quarkonium production, this is probably the less demanding one, as the formalism and the tech- 
niques to be used are a priori rather well established. Nevertheless there are some subtleties 
related to the regularization of the infrared (IR) and ultraviolet (UV) singularities which need to 
be addressed with some care. Furthermore, as pointed out in a preliminary account of part of this 
work the impact of NLO corrections can be significant and a general study of their effects is 
necessary. 

To carry out our calculations, we need a framework for calculating NLO inclusive production 
cross sections and inclusive annihilation decay rates. As well known, in perturbative calculations 
of the short distance coefficients beyond leading order in as, the most convenient method for 
regulating both UV and IR divergences is dimensional regularization. On the other hand most 
calculations of production cross sections and decay rates for heavy quarkonia have been performed 
using the covariant projection method which involves the projection of the QQ pair onto states 
with definite total angular momentum J, and which is specific to four dimensions. 

In this paper we generalize the method of covariant projection to D = A — 2e dimensions and 
perform all our calculations within this framework. Very recently a different method for calcu- 
lating production cross sections and decay rates, which fully exploits the NRQCD factorization 
framework and bypasses the need for projections - the so called "threshold expansion method" - 
has been generalized to D dimensions so that dimensional regularization can be used [^, We 
will show that, whenever common calculations exist, the results obtained using the two techniques 
coincide. 

The knowledge of perturbative corrections to quarkonium production is rather limited. The 
only examples of full NLO calculations we are aware of are the following: 



total hadro-production cross-sections for ^5*0 states 



inclusive px spectrum of J/ijj in photo-production, within the GSM 



corrections to the polarization of J/ip produced via gluon fragmentation |^ . 
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In this work we present the first calculation of the O(a^) total cross-sections for hadroproduc- 
tion of several QQ states of phenomenological relevance: ^S'o^'^', '^Si^"^ and '^Pj'^'^\ where the right 
upper index labels the colour configuration of the QQ pair. In addition, we repeat the calculation 
for ^5*0'^^ states within our formalism, in order to establish consistency with the previous calcula- 
tions of this quantity [BTl 0. We also calculate the O(a^) light-parton decay rates for ^Sq , ^S'l and 



^Pj colour-singlet and colour-octet states^ Some of these results are new, some have been known 
for a long time. In the case of ^Sq^^ decays we find agreement with previous results @, |^. In the 
case of ^P^j^ and ^S^q^ decays we find a disagreement with the previous calculations, reported in 
and in |^0|, respectively^. While these discrepancies have a negligible numerical impact, it would 



be interesting to see in the future new independent calculations performed, to confirm either of 
the results. 

This paper is structured as follows. In Section |^ we introduce our formalism, with emphasis on 
the extension to D dimension of the covariant projection method and on its link with the NRQCD 
factorization approach. Section ^ gives a brief general description of the NLO calculation. In par- 
ticular, we describe the technique used to identify the residues of the IR and collinear singularities 
and to allow their cancellation without the need for a complete D-dimensional calculation of the 
real-emission matrix elements. 

Section ^ describes the behaviour of the soft limit of the NLO real corrections, and Section ^ 
presents the results for the various components (real and virtual corrections) of the NLO result 
for the decay widths. Section ^ describes how, within the NRQCD formalism, the Coulomb 
singularity of the virtual corrections and the aforementioned infrared singularity which appears 
in NLO corrections to P-wave decays can be disposed of. Finally, Section ^ presents the NLO 
results for the production processes. 

Appendix ^ collects symbols and notations and Appendix ^ collects the results for the Born 
cross sections and decay rates in D dimensions. A summary of all results is provided in Ap- 
pendix |C|, where hadroproduction cross sections and decay widths into light hadrons are presented 
in their final form, after the cancellation of all singularities. In this Appendix we also comment 
on the discrepancies observed when comparing the results for ^Pq^', ci-nd ^Sq^' decays with 
those presented in previous calculations. Results for quarkonia photoproduction and decay into 
one photon plus light hadrons can be easily extracted from these calculations: explicit results will 
be presented in a forthcoming publication. 

Appendix presents, for comparison, a calculation of the NLO virtual corrections performed 
using the threshold expansion method. The results agree with those obtained using the D- 
dimensional projector technique. Appendix ^ presents the results for the g ^ Xj fragmentation 
function obtained within our formalism, showing again agreement with the results of the threshold 



expansion technique 42 



A preliminary account of some of the results contained in this paper, together with a first 
phenomenological study of NLO Xc b total production cross sections at fixed target and collider 



energies, was presented in ref. [36|. A more complete study, including the additional processes 



calculated in the present work, will be the subject of a separate publication. 



^Notice that in the case of colour-singlet '^Si and "^Pi states the 0{al) result is however LO only. 
^After this paper was released as a preprint, the authors of ref. |^ reviewed their calculation. We have been 
informed that their final result now coincides with ours. 
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2 Introduction to the formalism 



The calculation of cross sections and decay rates for heavy quarkonia deals with two kinds of 
contributions: short distance ones, related to the production or annihilation of the heavy-quark 
pair, and long distance ones, related to non-perturbative transition between the quarks and the 
observable quarkonium state. Rigorous theorems exist [|l5l proving the factorization of these two 



stages in the case of inclusive quarkonium decays. Less rigorous, but widely accepted, proofs have 
also been given of the factorization between the perturbative and non-perturbative phases in the 
hadropro duct ion case [^. In order to lay down the strategy of our calculation, and establish 
some notation, we briefly review here the essence of these factorization statements in the case 
of hadropro duct ion. Figure |I| shows a sketch of a typical diagram contributing to the inclusive 
production of a quarkonium state H . Factorization implies that we can identify an intermediate 
perturbative state composed of a pair of on-shell heavy quarks, whose non-perturbative evolution 
will lead to the formation of a physical hadron H plus possibly a set X of light partons. The set 
of light partons Y consists of additional states produced during the hard collision that created 
the QQ pair. Whether a light parton belongs to the set X or to the set Y depends on our 
(arbitrary) choice of factorization scale. The possibility to maintain the freedom to select this 
scale at all orders of perturbation theory is a consistency check of the factorization hypothesis, 
and leads to renormalization-group relations between the non-perturbative matrix elements which 
describe the long-distance physics |TB|. The factorization theorem states that the effect of soft 



gluons exchanged between the sets X and Y cancel out in the production rate, and their effect can 
therefore be neglected. The emission of hard gluons between X and y , furthermore, is suppressed 
by powers of the ratio between the scales of the soft and hard processes. 

At the amplitude level, and assuming the validity of the factorization theorem, the processes 
shown in fig. |1| can then be represented as follows: 

M = %j{QQ + Y) {HX\W,^^\0) (2) 

where Tij is the matrix element for the production of the QQ + Y final state with the heavy quark 
spinors removed, i and j are spinorial indices, and ip is the fermion field operator. Summing over 
a complete set of states containing the heavy quark pair, one then obtains: 

M = Y.iH^\QQM{QQ[^Mi^j\^)%j = Y.(Hx\QQ[n])Ti[rn^^^] , (3) 

n n 

where 77*^"^ selects the quantum numbers of the quark pair relative to the state n. Its specific form 
for S and P wave states will be given in the following. The first term on the right-hand side is 
interpreted as the overlap between the perturbative QQ state and the final hadronic state. It can 
be written as the transition matrix element of an operator , bilinear in the heavy-quark field, 
between the vacuum and the final state {HX\. We can therefore define the following quantity, 
following the notation introduced in p 



X 

where 



V"" = J2\HX){HX\ . (5) 



X 
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Figure 1: Diagrammatic representation of factorization in parton + parton H + X . 



We finally get the standard result WB: 



da{H + X) = Y,d^iQQ[n]+X){0"{n)) (6) 

n 

The quantity {0^{n)) is proportional to the inclusive transition probability of the perturbative 
state QQ[n] into the quarkonium state H. Notice that in principle we could have specified some 
detail of the inclusive state X, for example we could have specified the fraction z of light-cone 
momentum of the pair QQ[n] carried away by X, in which case we could have defined, starting from 
eq. (^), the equivalent of a non-perturbative fragmentation function (as opposed to an inclusive 
transition probability). Such a fragmentation function provides a more detailed description of 
the hadronization process, which can be used to parametrize potentially large higher-order f ^ 
corrections, such as those appearing near the boundary of phase-space. The factorization theorem 
allows us to extract it from some set of data, and its renormalization group properties should 
allow its use in different contexts. These more general non-perturbative quantities have also been 



introduced and discussed recently in ref. |3^, where several interesting applications have been 
explored. 

We now come to the discussion of the set of rules which define the perturbative part of eq. @, 
namely the projection over the perturbative states QQ[n\. These projections have been known 
for some time and have been used since then for the evaluation of quarkonium production 
and decays We will extend these rules here for use in dimensions D = 4 — 2e, in view of 
our applications to NLO corrections, where dimensional regularization techniques are most useful 
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to handle the appearance of infrared (IR) and ultraviolet (UV) divergences. With the exception 
of the normalization of spin and colour factors, discussed below, our choice of normalization is 
consistent with the standard non-relativistic normalizations used in []T^ to define the operators 



entering the definition of the non-perturbative matrix elements {O^ {n)). 

The spin projectors for outgoing heavy quarks momenta Q = P/2 + q and Q = P/2 — are 
given by: 

Ho = (^-^-"^iTsR + ^ + m), (7) 



m 




for spin-zero and spin-one states respectively. In these relations, P is the momentum of the 
quarkonium state, 2q is the relative momentum between the QQ pair, and m is the mass of the 
heavy quark Q. The above normalization of the projection operators corresponds to a relativistic 
normalization of the state projected out. Higher-order powers of g, not required for the study of 
S and P-wave states, have been neglected in eqs. and (||). 

The use of these operators, equal to those used in D = 4 dimensions, requires some justification. 
In D 4 dimensions, the product of two spinorial representations gives rise in fact to other 
representations, in addition to the spin-singlet (scalar) and spin-triplet (vector) ones. For example, 
the equivalent of a quarkonium state in D = 6 can have a total spin corresponding to a scalar, 
a vector, and a fully-antisymmetric three-index tensor of S0{5, 1). The total dimension of these 
representations is indeed 1+5-1-10=16. Once we move away from = 4, therefore, we should in 
principle take into account the existence of other states in addition to the S = and S = 1 ones 
found in D = 4. For non-integer D, in particular, we should deal with an infinite set of states, 
since the Clifford algebra becomes infinite-dimensional. The projectors defined in eqs. (|^,H) are 
still, nevertheless, the right operators needed to define the scalar and vector states. One can 
neglect the presence of higher-spin states in the infinite- mass limit, since spin- flip transitions, 
which would mix scalars and vectors to higher spins, vanish. With finite mass one can neglect the 
mixing among these operators provided one works at sufficiently low orders in f . At higher orders 
in V spin-flip transitions occur, and in principle mixing with higher-spin states should be handled. 
In the calculations illustrated in this paper this never occurs, and we can safely work with just 
the states we are interested in, namely the D-dimensional scalar and vector. 

The problem of dealing with the higher-spin states goes beyond the scope of this work, and will 
not be analyzed here in any detail. Nevertheless we want to provide some argument to suggest that 
even at high orders in v no problems are expected. The argument proceeds as follows: higher-spin 
operators should vanish in 4 dimensions, and therefore their effect should be suppressed by at least 
one power of e. In order to contribute to a cross-section in 4 dimensions, they must be accompanied 
by some 1/e pole. However IR and collinear poles cannot appear: IR poles arise from the emission 
of soft gluons, which cannot change the spin, and collinear poles admit a factorization in terms 
of lower-order amplitudes times universal splitting functions. Therefore no new operators can 
appear in the collinear limit, and higher-spin operators should decouple. As for UV poles, in the 
cross-section calculations presented here they are all reabsorbed by the standard renormalization 
procedure. We believe it should be possible to set the above arguments on a more solid footing. 
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We point out that a proof of the decouphng of higher-spin evanescent operators is also required 
when using the D-dimensional threshold expansion technique introduced in ref. [Q, since in its 



current formulation the explicit assumption is made that only spin-0 and spin-1 operators are 
relevant. 

The D-dimensional character of space-time is implicit in eqs. and appears explicitly 

when performing the sums over polarizations, as shown later. The manipulation of expressions 
involving 75 is carried out by using standard techniques, as discussed in the following. The colour 
singlet or octet state content of a given state will be projected out by contracting the amplitudes 
with the following operators : 

Ci = —7= for the singlet (9) 



Cg = V2T^j for the octet (10) 

The projection on a state with orbital angular momentum L is obtained by differentiating L 
times the spin-colour projected amplitude with respect to the momentum q of the heavy quark 
in the QQ rest frame, and then setting q to zero. We shall only deal with either L = or L = 1 
states, for which the amplitudes take the form: 

As=Q,L=Q = Tr [Ci7o^]|g^o ^P^^ singlet S states (11) 

As=i,L=o = CaTr [C ilf A]\g^Q Spin triplet S states (12) 

As=o,L=i = ^fi^r~'^^ P-^0'^]|g=o Spin singlet P states (13) 



As=i L=i = £a0^TT [C n^A]\ Spin triplet P states (14) 



A being the standard QCD amplitude for production (or decay) of the QQ pair, amputated of 
the heavy quark spinors. The amplitudes As^l will then have to be squared, summed over the 
final degrees of freedom and averaged over the initial ones. The selection of the appropriate total 
angular momentum quantum number is done by performing the proper polarization sum. We 
define: ^ ^ 

n„/3 = -^„/3 + -j-^ , (15) 

where M = 2m. The sum over polarizations for a state, which is still a vector even in 
D = A — 2e dimensions, is then given by: 

^e^e*, = n„a/ (16) 

In the case of ^Pj states, the three multiplets corresponding to J = 0, 1 and 2 correspond to a 

scalar, an antisymmetric tensor and a symmetric traceless tensor, respectively. We shall denote 

iJ) 



their polarization tensors by S^J . The sum over polarizations is then given by: 



1 

D 



^aP^S]^' - -UapUa'fS' (17) 
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Yl^ccp^a'P' = 2[naa'n/3/3' - Uap'Ua'p] (18) 



1 1 

'^^alS^a'P' = 2^""'"'^/^'^' ~'~ nQ/3'na//3] — jj _ ^a/B^a'p' (19) 



for the ^Pq, ^Pi and states respectively. Total contraction of the polarization tensors gives the 
number of polarization degrees of freedom in D dimensions. Therefore 

Nss, = E^ae; = n„„ = L> - 1 = 3 - 2e (20) 

for the ^5"! state and 

Nj-E^Ssi'J* (21) 

Jz 

for the ^Pj states, with 

iV.^ (^-^f-^) =(3-2.)(l-.). ;V.= (^ + ^f-^) = (5-2.)(l-.). (22) 

The apphcation of this set of rules produces the short-distance cross section coefficients a or 
the short-distance decay widths F for the processes: 

da{tj ^ 25+i^;i,8]) ^ ^El , (23) 

(25+i^;i,8l ^ ^ , (24) 

s being the partonic centre of mass energy squared and 2m representing the mass of the decaying 
QQ pair. To find the physical cross sections or decay rates for the observable quarkonium state H 
these short distance coefficients must be properly related to the NRQCD production or annihilation 
matrix elements (Off^^^i^^+'^Lj)) and {H\0[i,8]{'^^+'^Lj)\H) respectively. 

The above procedure contains some freedom in the choice of both the absolute and relative 
normalization of the NRQCD matrix elements. We could in fact decide to shift some overall 
normalization factor from the long-distance to the short-distance matrix elements. A natural 
choice for the absolute normalization of the non-perturbative matrix elements is obtained by 
requiring the short-distance coefficients to coincide with those which appear when the expectation 
value of the NRQCD operators is taken between free QQ states. The relative normalization 
between colour-singlet and colour-octet operators is obtained by imposing the decomposition: 

Otot - Ci + C»8 . (25) 

The identification of the colour-singlet and color-octet components is obtained by using the Fierz 
rearrangement: 

= + 2i;^^7^,., (26) 

The identification of the spin-0 and spin-1 components is obtained by using the Fierz rearrange- 
ment: 



(27) 
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with the normahzation of the Pauh matrices fixed to the canonical one: 



{a\ a''} = 26"'' . (28) 
Combining the two results, we obtain the following decomposition: 

V'V X'X = ^ (^^X XV + ^^cr*X X'^crV) 

+ ^p^T^'x X^T'^ip + ^Ij^T^a'x X^T^a'ijj . (29) 

We find it natural to use this decomposition to identify the normalization of the NRQCD oper- 
ators, which are presented in full detail for S and P-wave states in Appendix The resulting 
normalization for the colour-singlet part differs from the usual one found in the literature. Com- 
pared to the conventions of BBL |T5|, our normalizations are given byQ: 



iBBL 



O. ^ ^. (30) 



' c 

= O^"^. (31) 

We stress that, in addition to being more natural, (e.g. in the case of 5" waves the non-perturbative 
matrix elements for colour singlet states coincide exactly with the value of the total wave functions 
evaluated at the origin), our definition is probably more adequate when trying to estimate the 
order of magnitude of non-perturbative matrix elements using velocity-scaling rules. For example, 
it is known that the ratio between the matrix elements of the colour-octet and colour-singlet ^Si 
operators should scale like w^. The extraction of (Og '^^^('^S'l)) from the Tevatron data |2^, 0, pT] 
leads to values in the range 0.65 1.4 x 10~^ GeV~^. The value of {Of'^^^i^Si)) is given by 
1.3 GeV~^. The ratio between the two, of the order of 0.5 -h 1 x 10~^, is quite smaller than the 
estimated value of f ^ ~ 0.06. Using our normalization of the NRQCD operators, the ratio becomes 
3 ^ 6 X 10~^, which is consistent with the velocity scaling rule. This behaviour is confirmed by 
similar results obtained in the case of ip' and P-wave states. 

Having done this, the cross sections and the decay rates read: 

-(U - ^'-'l}'"' ^H) = ai^J -^^L;^'«l)i^M(!^M, (32) 

-^^ col-^^ pol 

T{H ^ ,^-) ^ tes+iL/'^^ ^ zj){H\Oi,,s]e'^'Lj)\H). (33) 

Ncoi and Npoi refer to the number of colours and polarization states of the QQ[^^^^Lj] pair 
produced. They are given by 1 for singlet states or A^"^ — 1 for octet states, and by the D- 
dimensional A^j's defined above. Dividing by these colour and polarization degrees of freedom in 

^Strictly speaking, consistency with the relativistic normahzation of the QQ state imposed by eqs. (|^) and (^) 
requires the NRQCD matrix elements to be evaluated using a relativistic normalization for the quarkonium state. 
The normalization factor 2M, implicitly assumed in the unitary sum over intermediate states used in eq. 
rescales however the matrix element to the value obtained using standard non-relativistic normalizations, up to 
corrections of order For this reason, in the following we will consider NRQCD matrix elements evaluated 

using non-relativistic normalizations, as customary, and no additional 2M normalization factor is required. 
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the cross sections is necessary as we had summed over them in the evaluation of the short distance 
coefficient a. Analogously, these degrees of freedom are summed over in the decay matrix element, 
but averaged as part of the definition of the short distance width F, when mediating over the QQ 
initial state. 

The matrix elements appearing in the equations above are of course meant to be the bare 
D-dimensional ones whenever D-dimensional cross sections or decay rates are to be calculated. 
Making use of their correct mass dimension, 3 — 2e and 5 — 2e for S and P wave states respectively 
(see Section |^), gives the right dimensionality to /^-dimensional cross sections and widths, i.e. 
2 — D = —2 + 2e and 1, respectively. 



2.1 Digression on the use of 75 in dimensional regularization 

Once the right counting of degrees of freedom is taken into account via the polarization sums 
eqs. (|T^,0,|T^) , further potential inconsistencies which could, a priori, spoil the possibility of 
using projectors on D-dimensional amplitudes, are related to the the definition of the matrix 75 
in arbitrary dimensions. Several prescriptions have been introduced and shown to be consistent 
in specific problems. The most common ones are the naive dimensional regularization (NDR) 
and the t'Hooft-Veltman dimensional regularization (HVDR) p4 |. 

Within the ffist prescription, 7^ is defined by the property that it anti-commutes with all 7^^ 
in D dimensions 

{75^,7^ = 0, /i = 0,l,---,Z}-l, (34) 

and it satisfies 

(75^)' = 1 ■ (35) 



while, in contrast to (|3^), in the t'Hooft-Veltman prescription it holds 

^ nVi'i' , (36) 

and 

{75'''', 7^ = 0, /i = 0,l,2,3, (37) 

75^^,7'^] = 0, /x = 4,---,D-l. (38) 

Since it is defined explicitly by construction in eq. (p6D, •y^^ is unique and well-defined, while 
to this respect, 7^, is an ambiguous object. For example, while the HVDR scheme allows to 
recover the standard Adler-Bell-Jackiw anomaly, the NDR does not without additional ad hoc 
prescriptions. 

We do not expect any ambiguity to arise in our case, since in QCD parity is conserved and 
since no 75 insertions are present inside the loop amplitudes we will be evaluating. Extra care 
should instead be applied when dealing with, for example, corrections to the decay b Qcs s. To 
verify the independence from the scheme, we did the calculation of virtual diagrams contributing 
to NLO decay widths and cross sections of ^S"!)^'*^' states in both in the HVDR and NDR schemes. 
The results were found to coincide exactly, diagram by diagram. 
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Finally we would like to compare our technique with the ones available in the literature. 
Braaten and Chen determine the short- distance coefficients via the matching technique, 

generalizing to D — 1 spatial dimensions their threshold expansion method. This technique is 
briefly presented in Appendix and has been used in this paper as a check for the virtual 
contributions in the NLO calculations. Within this approach a quantity that is closely related 
to the cross section for the production of a QQ pair with total momentum P is calculated using 
perturbation theory in full QCD and expanded in powers of the relative 3-momentum q of the 
QQ pair and in a few non relativistic matrix elements. To exploit this procedure one is forced to 
choose a representation for 7-matrices. Braaten and Chen write 



f = ' . n , t = l,...D~l. (39) 





Once this choice has been made, one easily realizes that no room is left for selecting a 75 different 
from 

- f ; J ) . (40) 



which, as a simple calculation shows, satisfies eq. ([5^). This observation shows that in the D- 
dimensional threshold expansion method an implicit choice for dealing with 75 is made and it 
exactly corresponds to 7^*" = 7^^. 



3 General description of the calculation of higher-order 
corrections 

In this section we briefly describe the strategy of the calculation of higher-order corrections to 
decay widths and to total production cross sections. The reader who is not interested in the details 
of the calculation, can just read this section to get an idea of the general framework, and can skip 
the following sections where all the required calculations are carried out. 

A consistent calculation of higher-order corrections entails the evaluation of the real and vir- 
tual emission diagrams, carried out in D dimensions. The UV divergences present in the virtual 
diagrams are removed by the standard renormalization. The IR divergences appearing after the 
integration over the phase space of the emitted parton are cancelled by similar divergences present 



in the virtual corrections, or by higher order corrections to the long-distance matrix elements . 
CoUinear divergences, finally, are either cancelled by similar divergences in the virtual corrections 
or, in the case of production, by factorization into the NLO parton densities. The evaluation of the 
real emission matrix elements in D dimensions is usually particularly complex, and is presumably 
the main reason that has prevented the so far the calculation of NLO corrections to the production 
of P-wave states. In this paper we follow an approach already employed in |^ , whereby the struc- 



ture of soft and collinear singularities in D dimensions is extracted by using universal factorization 
properties of the amplitudes. Thanks to these factorization properties, that will be discussed in 
detail in the following section, the residues of all IR and collinear poles in D dimensions can be 
obtained without an explicit calculation of the full D-dimensional real matrix elements. They 
only require, in general, knowledge of the D-dimensional Born-level amplitudes, a much simpler 



12 



sJULJULSUr 



OLfiJLSUU 



OLOJLfiJ 




D2 D3 

Figure 2: Diagrams for tlie qq and gg Born amplitudes. 



task. The isolation of these residues allows to carry out the complete cancellations of the relative 
poles in D dimensions, leaving residual finite expressions which can then be evaluated exactly 
directly in D = 4 dimensions. In this way one can avoid the calculation of the full D-dimensional 
real-emission matrix elements. Furthermore, the four- dimensional real matrix elements that will 
be required have been known in the literature for quite some time 
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4 Soft emission behaviour 

We discuss in this section the factorization properties of the real-emission amplitudes in the soft- 
emission limit. The factorization formulae presented here will be used in the following sections to 
isolate the IR poles and cancel them against the singularities of the virtual processes. 



4.1 Soft factorization in 3g amplitudes 

We start by considering the case of decays into gluons. At the Born level, the relevant diagrams 
are shown in fig. ^. The decay amplitude (before projection on a specific quarkonium state) can 
be written as follows: 

^Born = {ab)ij (Di + D3) + (ba),, {D2 - Ds) , (41) 
where we introduced the short-hand notation: 

(a... 6),, = (T,...T,),. , (42) 



13 



where the T" matrices are normahzed according to Tr(T'^T^) = 5°"^ /2 (more details on our conven- 
tions for the colour algebra can be found in Appendix |^. The terms Di, D2 and correspond 
to the three diagrams appearing in fig. ^ with the colour coefficients removed. Using this notation, 
the amplitude for emission of a soft gluon with momentum k and colour label c can be written as 
follows 



A. 



soft 



= g{cab)ij 
+ g{ahc)ij 
+ g{hca)ij 



Qk 

bee Qe 
bk~^ 



aec 
ak 



Qk 



atr 



bk ak 



(Di + D3) + g{acb),j 
{Di + Ds)+ g{cba)^J 
{D2 - D-i) +g{bac)ij 



bee 
bk 

bee 
Qk bk 

aec Q(^c 



aec 
ak 

Qec 



ak Qk 



{D2 - D^) 



(43) 



where Q and Q are the momenta of the heavy quarks, and a, b indicate the momenta and colour 
labels of the final state gluons. 

In the case of S'-wave decay we can set Q = Q = P/2. All terms proportional to Qec or 
to Qec vanish in the transverse gauge defined by the gauge vector {ko,—k). The soft-emission 
amplitude becomes: 

ber 



A 



soft 



9 



+ 9 



aec 



([«>c]6)ij — - {a[c,b])ij — 



{[b,c]a] 



bee 
bk 



{b[c, a]) 



aec 
ak 



{D2 - D,) 



(44) 



In the case of ^P-wave decays we also need the derivative of the decay amplitude with respect to 
the relative momentum of the quark and antiquark. In the soft-gluon limit, we obtain: 



dA, 



soft 



dqo 



q=0 



9 



+ 9 



i[b,c]a)ij^- {b[c, 



+ 2gj^{[iabc),, + icab),,]iDi 



d{Di + D^) 
dqa 

d{D2 - Ds) 
dqa 

Ds) + [ibac)i, + {cba)ij] {D2 - D3)} (45) 



It is easy to prove with an explicit calculation that the two terms proportional to e" give vanishing 
contributions, in the soft-gluon limit, when projected on all ^Pj states. Therefore the structure 
of the soft limit of the derivative of the amplitude is similar to that of the amplitude itself. 

We now proceed to prove explicitly the factorization of the soft-gluon emission probability. 
We will indicate with {i = 1, 2, 3) the projection of the diagrams on a given quarkonium 
state, and with ABom and Agoft the projected amplitudes for the Born and soft-emission processes, 
respectively. 

The colour coefficients for the decay of a colour-singlet state can be obtained by using the 
projector SySki/Nc. The colour-singlet amplitudes therefore become: 

1 



A 



1] 

Born 



2N, 



■6abSij{Di + D2) 



(46) 
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A. 



soft 



aec bee 
ak bk 



{D1 + D2) 



and simple colour algebra gives the final result for the soft-emission: 

2ab 



col, pel 



[1] 12 
soft I 



The colour-octet case can be obtained from the following relation: 

col,pol col,pol col,pol 



(47) 



(4J 



(49) 



The first term on the right-hand side can be easily obtained with standard colour algebra, and 
the final result is: 



/ . I ^soft I 



coljpol 



2ab 



2ab 



{ak){bk) {Pkf 



/ . I^Bornl 



(50) 



col,pol 



4.2 Soft factorization in light-quark processes 

An approach similar to that used in the previous subsection can be applied to the case of soft- 
gluon emission in processes where the QQ pair Q decays into light-quark pairs (Q — > qq). The 
Born-level amplitude (which is non-vanishing only in the colour-octet case) is given by: 



A 



Born 



Do 



(51) 



where {i,j) and {k, I) are the color indices of heavy and light quark pairs, respectively, and where 
Do is shown in fig. |^, with the colour factors removed and before projection onto the quarkonium 
quantum numbers. The amplitude for emission of a soft gluon with momentum k and colour label 
c can be written as follows: 



A, 



soft 



9D0 
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+ -f5.k 


q^c 
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qk 


qk 


Qk 



Tl 
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kl 
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+ 2iv/- 
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1 



The colour-singlet and colour-octet projections are then given by: 



A 



soft 
[8] 



9 Do rpc. 



Qec Qec 



Qk Qk 



soft 



^soft - A. 



[1] 

soft 



(52) 

(53) 
(54) 



We study first the case of S-wave decays. As before, we can set to zero all terms proportional 
to Qec and Qec by choosing a transverse gauge. The colour- singlet amplitude obviously vanishes, 
and the colour-octet one reduces itself to: 



A, 



gDo 



soft 



Nc 



qec 
qk 



+ 



Ttj^ik 



1 



(55) 
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An explicit calculation gives the following factorization formula, which is non-vanishing only in 
the ^Si case since the ^Sq Born amplitude vanishes: 



soft I 



col,pol 



y I^Bornl 
col,pol 



2qq 



iqk){qk) 



2M2 Ca 



(56) 



where M = 2m. 



In the case of P waves we need to consider as well the derivative with respect to the relative 
momentum of the heavy quarks: 



'dA. 



soft 



dqa 



g=0 



(57) 



The terms proportional to dD^/dqa vanish because Dq corresponds to an S-wave process. Denoting 
by C the overall colour coefficient in the previous equation, choosing a transverse gauge where 
ec ■ P = and using the projection operators for ^Pj-waves given in section ^ we obtain the 
following result for the quarkonium-decay amplitude: 



A. 



soft I p. J 



fl'-^C ABorn(^5'i(eefr)) • 



(5J 



The amplitude was written in terms of the amplitude for the production of a ^5*1 colour-octet 
state, with an effective polarization eefr given in terms of the polarizations of the soft gluon and 
of the ^Pj state as follows: 

(59) 



-eff,/3 



The correlations between the polarization of the gluon and of the Q state induce a dependence 
of the soft-emission amplitude on the relative direction of the gluon and of the light quarks. 
Nevertheless, if we average over the soft gluon D — 1 spatial directions we get: 



[d^^-^ D-2 

pol 



(60) 



and one can easily compute the sum over effective polarizations 



fdnr'^j J* _ yv^-^n (P) 



D-1 



(61) 



Here Nj is the number of degrees of freedom of the ^Pj state in D dimensions. With the above 
expressions at hand, it is straightforward to square the amplitude in eq. (EW) and obtain 



dn^-'Y.\A,oitCPj)\' = n 



D-iiD-2) 



colours 



(Pk)' 



iVjEl^BornC^l' 



3 q[S]\\2 



(62) 



The above formula embodies the well known singular behaviour of the P-wave decay into light 
quarks [|10|, to be absorbed into the colour-octet ^Si NRQCD matrix element, as shown in ref. 
and discussed in Section 0. This formula also shows that this phenomenon is present not only for 
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colour-singlet P wave states, but for colour octet ones as well. The colour factors (averaged over 
initial states) are given by: 



E IC'^'P = Cf, (63) 



colours 



E l^'^'r = Bp, (64) 



colours 



for the colour- singlet and colour-octet cases, respectively. Cp and Bp are defined in Appendix 0. 

5 Decay Processes at NLO 

The previous study of the soft behaviour of the real-emission amplitudes will now be used to 
calculate the full set of NLO corrections to decay and production of quarkonium states. The 
general factorization formulae that we proved in the previous section will be useful to efficiently 
extract the singular behaviour of the real-emission amplitudes, using the general formalism of 



refs. |^5[. In conjunction with the universal behaviour of amplitudes in the coUinear regions, the 
knowledge of the residues of the soft singularities will enable us to extract all poles of infrared 
and collinear origin, cancel them against the poles in the virtual amplitudes and in the parton 
densities, and be left with finite terms which only depend on the real-emission matrix elements 
in 4 dimensions. This technique saves us from the complex evaluation of the full matrix elements 
for real emission in D dimensions. In this section we confine ourselves to the case of decay rates. 
We will discuss the production cross-sections in Section |^. 

5.1 Real emission corrections in D dimensions 

5.1.1 Kinematics and factorization of soft and collinear singularities 

Let us consider the three-body decay processes Qt^'^] — > ki + k2 + k3, where Qt^'^] = QQ[^^^^Lj^^'^\ 
They are described in terms of the invariants 

^i = ~J^^ E^i = 2, (65) 

where P is the momentum of the decaying QQ pair of mass M = 2m. The three-body phase 
space in D = A — 2e dimensions is given by: 



N ^ 

^{2)jtII(^ - ^^)~' dx,5{2 -Y^x,) , (66) 

i=l 



where $(2) is the total two-body phase space in D dimensions: 

1 fAnV T(l-e 



6TT 
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and and K are defined as 

= r(l + e)r(l-e) ~ l + eV • (69) 

o 

For future reference, it is useful to define as well tlie following function: 

= {jfj ni + e), (70) 

where fi is the renormalization scale. 

Collinear and soft singularities in the matrix element squared 

M{Xi,X2,X3) = ^|A(xi,X2,X3)P (71) 

appear as single poles in either of the terms 1 — Xj. As a consequence, the function 

3 

fiXi,X2,X3) = YI{1 - Xi)M{xi,X2,X3) (72) 
i=l 

will be finite throughout the phase space. In terms of the function /, the differential decay width 
can be written as: 

^(2) N 



dT = ^ n (1 - ^«)"'"' dx, 5(2 - J2 fi^i, ^2, xs) . (73) 



i=l 

It is useful to introduce the variables (x, y, z), defined by: 

Xi = 1 — xy (74) 
X2 = I — x{l — y) (75) 
X2, = z (76) 



With these variables 



n(l-x,) = x\l-x)y{l-y) , (77) 

i=l 
3 

W_ dxi5{2 — '^Xi) = xdxdy, (78) 
1=1 

with z = X. The total decay width can then be written as: 

r = ^^y^'d^l'dy^~'~"'(^-^y'"iy(^-yr'"fix,y) (79) 

where 5 is a symmetry factor, equal to 3! in the case of 3-gluon decays, and equal to 1 in the case 
of gqq decays. The function 

f{x, y) = f{xi = 1 - xy, X2 = 1 - X + xy, X3 = x) (80) 

is finite for all values of x and y within the integration domain. Therefore all singularities of the 
total decay rate can be easily extracted by isolating the e poles from the factors in eq. ([79|) 
explicitly depending on x and y, as discussed in the next subsection. 
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5.1.2 Colour-singlet three-gluon decays 



For the sake of definiteness, we present now the detailed calculations relative to the decay of a 
colour singlet state into three gluons. The colour octet can be worked out in a similar manner, 
and will be given with fewer details in the following. 

To exploit the total symmetry of the three gluons in the final state we restrict the integration 
to the domain < xi,X2- In terms of the variables x and y, this corresponds to 



2 1 1 
0<x<- , 2--<y<--l. 

3 XX 



(81) 



This choice corresponds to integrating over one third of the entire phase space, multiplying the 
result by the multiplicity factor 3. To proceed, we make use of the following expansions vahd for 
small e: 

-1-26 _ , ^1^ ^^/^logx' 



X 



4'' /I 



1/2 V ^ J 1/2 



where we introduced the "1/2" distributions defined as follows: 

f.l/2 1-1/2 



dx [d{x)]^l^t{x) = ^ dxd{x) [t{x)-t{Q)] 



We then obtain: 
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log X 
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1/2 J 



(82) 



(83) 



+ 0{e' 



B4) 



The total decay width can then be written as the sum of three terms: 

r = ^x=G + ro<a;<l/2 + ^l/2<x<2/Z , 



where 
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l/2<a;<2/3 
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2M K S Ji/2 



/•■Z/6 rl/x — l 
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/x x{l-x)y{l-y) 



(85) 



(86) 



(87) 
(88) 



Notice that the last term is finite, since no pole is present within the integration domain. It can 
therefore be computed directly in 4-dimensions. 
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We shall now proceed to the evaluation of each of these terms, starting from the first one. The 
X — s> limit is the soft limit for k^. In this limit, eg . (^81) gives: 



Mix,y) 



M 



Born 



9) 



M2 x^y{l - y) 

where af^ = a^/i^^ is the dimensional coupling constant in D dimensions and A^Born is the Born 
amplitude squared and averaged in D dimensions for the Q ^ gg processes. As a result, 

4^ 1 $(2)A^ 
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dy [yil-y)] ^ ' f{x = 0,y) 



2e2M 2! K 
1 $(2) levrCAtt^'') 



2M 2! K 
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(90) 
(91) 



In the cases in which A^Bom = 0, such as for example Q = ^Pi, AlBom = and the left-hand side 
of the above equation, together with r^_o, vanish. 

In order to evaluate ro<a;<i/2, we use the following relation: 



-l-e 



1 



[S{y)+Sil-y)] + 



1 



as follows: 



+ Oie) 



(92) 



With obvious notation, we can therefore decompose ro<a;<i/2 

^0<x<l/2 = ^y=0 + ^y=l + Tfinite • (93) 

The limits y ^ and y ^ 1 correspond to the coUinear limits /c2 || ^3 and ki \\ k^, respectively. 
In these limits we can apply the universal factorization of coUinear singularities 



f{x,y = 0) = f{x,y = l) 



Pgg{x) MBovnX{l - X) 



where Pgg{x) is the D-dimensional Altarelli-Parisi splitting kernel: 



X \ - X ,^ , 
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— X X 
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and obtain the following result: 
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We can now collect all the results obtained so far in the following expression 



where 
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^Born/e(M^) (- + - + - --7r^ + - log 2 



(9J 
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is a universal term, whose dependence on the quarkonium type is only implicit in the factor FBom, 
and where 



r 



finite 



1 $(2) A^ 

2M 2! K 



dx 



dyf{x,y) 




+ ri/2<..<2/3 (99) 



is a finite expression, explicitly dependent on the quarkonium state through the specific form of 
the function f{x, y). In particular, this is the only term which survives in cases where A^Bom = 0, 
consistently with the form taken by eq.(p9|) in the limit e — 0. Since Ffinite is free of soft and 
coUinear singularities, it can be evaluated directly in four dimensions. In particular, this implies 
that the matrix elements for the real process Q ggg can be calculated in D = 4 dimensions, 
with a significant reduction in the complexity of the calculation. The evaluation of the Ffinito 
contributions is a tedious but straightforward calculation, which makes use of the 4-dimensional 
matrix elements given in ref. for the colour-singlet states. The explicit evaluation leads to the 
following result (valid for FBom 7^ 0): 
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(101) 
(102) 
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The processes ^5*1 ggg and ^Pi — > ggg are completely finite. In absence of real or virtual 
contributions from two-gluon final states, they provide the LO contribution to the gluonic decay 
rates: 



999) 
999) 
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^ ^ TT^ {H\0,{'S,)\H) 
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(104) 
(105) 



These results agree with previous evaluations (see e.g. Schuler [0). 

We would like to point out an interesting fact related to the size of the non-universal finite 
corrections appearing in eqs. ([Toll) through (|10|): although the size of the contributions from the 
rational number and from the term proportional to vr^ are large in each case, there is always a 
cancellation between them which reduces their sum to about 1/10 of their individual value. The 
same accurate cancellations take place in the processes that will be considered in the following 

fl 8l 

subsections. In all cases the finite coefficient An is a number of order 1. 



21 



5.1.3 Colour-octet three-gluon decays 

The case of colour-octet decays can be analysed in a similar way, the only difference being in the 
soft limit of the amplitude squared, given by eq.(|50|): 



M{x,y) — A^Bo.n, (106) 

where A^Bom is the Born amplitude squared and averaged in D dimensions for the Q^^l gg 
processes. Following the same steps as above, we get: 

r(Q[8]^3^?) = rSi^ + Tfinite, (107) 

where 



r'l = ^^re„„/.(Af^) + 1 + f - + lllog 2) . (108) 



is the universal term for colour octet decays. Once the finite part is calculated, making use of the 
4-dimensional amplitudes derived in ref. 1^^, we obtain the following results: 



r(Q[«l - ggg) = ^T^o.M'^ - 99)UM') + l + A^) , (109) 



with 

1 n/i 

TT^ , (110) 

a':l = ^-^7i\ (111) 

Afl = ni:i-—Tt^. (112) 

^2 432 72 ^ ' 

The processes '^Pi^*^ — > ggg and ^Si-^^ ggg are completely finite because of the vanishing of 
two-gluon amplitudes. They read: 

r-CPr ^ ^ a^. (^ - f .') <MM> , („3) 
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5.1.4 Colour singlet qqg decays 

The ^Sq^^ — » qqg decay rate can be calculated along the same lines as above, starting from eq. (|79D . 
No soft-gluon pole is present for x — 0, and the collinear singularities are extracted by using the 
Altarelli-Parisi factorization. The result is given by: 
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The process ^Si^'^ qqg is equal to zero. The resuhs for P-wave colour-singlet states are given 
by: 



where 



-g''^^^"^ 
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and /e is defined in Appendix 0. 

In the above equations the 1/e poles proportional to Tp are due to the gluon splitting into 
coUinear quarks. They will be cancelled by the virtual corrections to the two-gluon decay process. 
The singular J-independent pieces, proportional to C^, come from the soft-gluon pole isolated in 
eq. (0), and will be cancelled by the addition of the ^^i'^^ qq contribution to the decay width 
(see Section ||). 



5.1.5 Colour octet qqg decays 

The process ^^o^^^ qqg is completely analogous to the colour singlet one: 
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The process Si 



3 c [8] 



qqg shows collinear and IR singularities that are expected to cancel 



against the virtual ones: 

T^C'Si'^'^q-qg) = rf_(35j^Ugg)%,(M^){Cp(l + A 



Ca— 
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(119) 



The colour-octet P-wave decays exhibit the same singularity structure as the colour-singlet 
case, and are given by: 



r^(3p|^l qqg) = n,rf_(3p] 
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and the quantities aj are defined in ( |117| ). Bp is defined in Appendix |^. 
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5.2 Virtual corrections 



We present in this section the results of the calculation of the 1-loop diagrams necessary for the 
evaluation of the virtual corrections to the production and decay matrix elements. The final 
results for the decay rates, obtained by combining the results of this section with those on real 
emission presented in the previous sections, are collected in Appendix |C[ 

The explicit calculation of the virtual contribution Fy was done using dimensional regulariza- 
tion, both for the UV and the IR divergences. We have introduced the parameter v defined in the 
Appendix ^ in order to regularize the Coulomb singularity. The quantity v represents the velocity 
of the heavy quarks in the quarkonium rest frame, and is kept small but finite. The Coulomb 
singularities appear as poles in the relative velocity, i.e. as l/2v. The relevant Feynman diagrams 
are shown in figures ^ and H, and the results are given diagram by diagram in tables |l|, 0, |^ and 
|[ In these tables we report the contribution of each diagram fc, indicating separately the colour 
factors fk relative to the colour singlet and colour octet cases. 

The expressions Vk appearing in the tables are defined by the following equation: 

Tv = rBorn-MM')Y.V,fk , (121) 

^ k 

where the sum extends over the set of diagrams. 

The virtual corrections are both infrared and ultraviolet divergent. While the infrared poles 
are cancelled when adding the real corrections, the ultraviolet ones are instead disposed of via 
coupling constant and heavy quark mass renormalization. The latter has already been explicitly 
performed in the on-shell scheme, which amounts to replacing the bare mass by 



—Cf 1E + log47r + - 
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The coupling constant renormalization has instead been left undone, to give in principle anyone 
the possibility to perform it in a scheme of his own choice. We will do it in the MS scheme through 
the replacement of the bare coupling constant by 



(— 7£; + log47r) 

2n Ve,Tv / 



(123) 



We point out that we carried out the evaluation of the virtual corrections using two independent 
techniques. In addition to using the D-dimensional projection operators constructed in sec. |[ 
we also used the threshold-expansion technique introduced in ref. [0, |2|, and reviewed here 



in Appendix 0. The results obtained in the two cases coincide exactly, diagram by diagram, 
providing an important consistency check of our calculation showing the equivalence of our D- 
dimensional projections to the threshold expansion method. 

5.2.1 Colour-singlet gg virtual decays 

The singularity structure of the virtual corrections is dictated by the renormalization properties 
of the theory, by the universal form of the Coulomb limit, and by requirement that soft and 
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Figure 3: Feynman diagrams contributing to the one-loop corrections to the processes 
gg[3p[i.8l] ^ gg^ QQppj^'^]] ^ gg and QQ[i4'''l] - gg. 



coUinear singularities cancel against the real corrections evaluated above. The form of the virtual 
corrections to the decay rate is therefore the following: 

r?(S™ -gg) = rg„„(Ql'l - 99)^/,(Af^) 

where we explicitly labelled the e's to indicate their origin, and where all of the state dependence 
is included in the finite factor Bq. The heavy quark-antiquark relative velocity 2f , fog and /^{M"^) 
are defined in Appendix 0. Uf is the number of flavours lighter than the heavy, bound one, as a 
consequence of the heavy quark never appearing in the virtual loops. 

Summing the contribution of all diagrams, we obtain the following results for the colour singlet 
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Table 1: First column: diagram labels according to the conventions of fig. ^ Second 
column: diagram-by-diagram contribution to the virtual corrections to the matrix element 
squared for the processes QQ[^5q^'^'] — > gg and with colour factors set to 1. Third column: 
colour factors relative to the gg colour singlet decay. Fourth column: colour factors 
relative to the gg colour octet decay. 



coefficients i?n 



< = Cf(-5 + ^)+Ca(i + ^,^), (125) 



7 7r2\ ^ fl 5 

L2 

1 5. TT^^ 



Bfl = CF{-i; + ^]+CA{^ + i^n'] , (126) 



3 4 ; V3 12 

'1 5, TT^ 

We recall that the virtual corrections to processes which are absent at the Born level vanish 



= -4CF + CA(^ + ^log2 + — ) . (127) 
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Table 2: Same as table |l], for '^P^ ' states. 



5.2.2 Colour-octet virtual gg decays 

Summing the contributions in tables || and 0, we get the following result for the colour-octet 
gg decays: 



X 



71 





where the state-dependent finite parts are given by: 

< = Cf (-5 + ^) + (2 + Ivr^) , (129) 
<^C^.(44).C^.(S + ilo.2.1.^), (130) 

< = -4C^ + (1 + ^ l°g2 + K^) . (131) 
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Table 3: Same as table ^ for '^P^' ^ states. 



5.2.3 Colour-octet qq virtual decays 

The only relevant process in this channel is relative to the '^S'f configuration, and the diagrams 
are shown in fig. ^. The contributions of the various diagrams are given in table ^ and their sum 
is 

TT 

with 

BPM^I] = Cp (-8 + ^vr^) + (f + ^log2 - ^) - '-^n,Tp . (133) 

6 Cancellation of IR singularities within NRQCD 

Throughout this paper, we witness the cancellation or removal of both infrared/collinear or ultravi- 
olet singularities by the usual QCD mechanisms: soft/virtual cancellation, coUinear factorization. 
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Figure 4: Feynman diagrams contributing to the one-loop corrections to the process 



mass and coupling constant renormalization. Two exceptions exist: the 1/v Coulomb singularity 
in the virtual corrections, see eqs. ( |124| , |128| , |132| ), and the 1/e infrared singularity which appears 



in the ^pj^'^l qqg ^eal correction processes in eqs. (|116| , |120|) (or in qq — pf'^^g, as discussed in 
the next section) cannot be eliminated via these standard mechanisms. Their removal is strictly 
related with the NRQCD factorization approach to quarkonium production and decay. It is within 
this approach that one finds a rigorous solution to this problem, previously dealt with in an empiric 
way, by absorbing the Coulomb singularity into the Bethe-Salpeter wave function and cutting off 
the infrared singularity with the binding energy of the quarkonium. 

Within NRQCD, one can determine the short distance coefficients of the various operators 
by performing a matching between cross sections calculated in perturbative QCD and perturba- 
tive NRQCD. Since, by definition, the two theories are equivalent in the long distance regime, 
all singularities of infrared origin appear equally in both calculations, and hence cancel in the 
matching. 
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Table 4: Diagram- by-diagram contribution to the virtual corrections to the matrix 
element squared for the processes QQ^S^\ —>■ qq. Diagram labels according to the 
conventions of fig. |^. 

Our approach does not perform an explicit matching, but makes use of projection tech- 
niques and then complements the short distance cross sections and decay widths with the non- 
perturbative NRQCD matrix element describing the transition to the observable quarkonium state. 
Rather than seeing the Coulomb/infrared singularities disappear in the matching, we will therefore 
cancel them by taking into account radiative corrections to the NRQCD matrix elements. 

Let us first describe with some details the cancellation of the Coulomb singularity. We discuss 
for definiteness the case of decays, but the same applies to production cross sections. The Coulomb- 
singular part of the virtual correction to a two-partons decay of a Q = QQ[^^^^Lj] state reads: 

ty{Q ^J) = tn^^^^f^(M') (c^ + ..) . (134) 

71 \ ZV J 

The coefficient C is a colour factor which, for Q in a colour-singlet or octet state, takes the value 

= Cp, 

C[8] = Cp-^Ca. (135) 

The above expression ( |134D , multiplied by the NRQCD matrix element {H\0{Q)\H) , provides 
the final answer for the virtual corrections to the decay width of the physical quarkonium state H. 
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Figure 5: Diagrams for NLO corrections to the NRQCD matrix ele- 
ments. 



However, to get rid of the Coulomb singularity, we must first evaluate in NRQCD the radiative 
corrections to this matrix element. This was done for the first time in ref. [0, and we reproduce 
their argument for illustration. 

{H\0{Q)\H) is depicted in fig. ^(a) in lowest order, and the higher order corrections important 
for the Coulomb singularities are shown in figs. ^(b,c). Let us consider, for instance, diagram ^(b). 

Working in Coulomb gauge, the gluon propagator is given by the sum of a transverse part. 



G transi}^ 



'1-5,q){1-5, 



and of a Coulomb part. 



(136) 



(137) 



It is this second part responsible for the singularity we are looking for. 

Inserting this term only and making use of NRQCD Feynman rules we get for diagram ^(b) 

d^k I I 1 



2m 



le 



(138) 



with the on-shell condition po = p^/2m, p being the heavy quark momentum and k the gluon one. 
The integral over k^ can be performed by contour integration, closing the contour in the lower 
half of the complex k^ plane and picking up the pole at fco = — Po + (p + k)^/2m — ie. The result 
is 

'^ = ^-"^7 (2^k^ k^ + 2p.k-.e - 

Being divergent, this integral must be performed within some regularization technique. Using 
dimensional regularization, and defining v = p/m, v = |v|, we find 



Av 



[1 + infrared singular imaginary part) 



(140) 
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From the diagram ^(c) we find Jc, similar to h but with opposite imaginary part. Summing 
the two contributions, the overall Coulomb correction to the matrix element is therefore: 



{H\0{Q)\H) = {H\OmH)Born [l + C 



.TTOs 

2v 



(141) 



The colour factor C evaluates to the same defined above in ( |135| ), according to the operator we 
are correcting (the central blob in figure |^) being a colour singlet or a colour octet one. 

For the decay width we have then 



rBorn + Tv{Q^ ij) + 
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Upon inspection, we can see that the singular 1/v terms cancel. Since no finite parts are introduced 
in this process, we can just drop the singular Coulomb term whenever it appears in the virtual 
contributions to cross sections and decay widths. 

Next we consider the problem of the infrared singularity which appears for instance in the 



3pji.8] _^ corrcctiou process. We have mentioned in Section [4.2| how the residue of this 

singularity is proportional to the matrix element squared for the process ^Sf"^ — *• qq, making it 
possible to absorb the singularity in the colour-octet wave function which accompanies this process. 
Again, NRQCD puts this cancellation on rigorous grounds. Radiative corrections to the colour 
octet matrix element {H\Os{^Si)\H) display a singularity which matches the one appearing in the 
P-state decay calculation. When summing P-wave and ^5*1 octet processes, the singularity cancels 
exactly. Some finite parts are however introduced by this process, which makes it necessary to 
carefully evaluate the diagrams involved. 

Let us then consider the radiative corrections depicted in diagrams |](d,e,f,g). It is this time 
the transverse part of the gluon propagator to be relevant, and the loop integral reads 
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Contour- integrating over ko around the ko = \h\ — ie pole we find 



d^k p-p'- (p-k)(p'-k)/k2 
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This integral is infrared divergent but ultraviolet finite. However, it has been argued [^, ^ 
that the correct way to perform it is to first expand [|15| the integrand in powers of p/m, p'/m 
and k/m, since it is for small values of these momenta that NRQCD is valid (though alternative 
approaches exist, see for instance ref. pHI). Such an expansion produces 
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which is now both infrared and ultraviolet divergent. We can perform it by using dimensional 
regularization and obtain 

where and euv are poles of infrared and ultraviolet origin respectively. It is to be noted 
that this ultraviolet divergence arises within NRQCD, and has nothing to do with the ultraviolet 
divergences which appear in the virtual corrections within full QCD, to be subsequently removed 
by renormalization of the coupling constant. 

The other three diagrams give identical result but with different colour factors when including 
the colour structure of the Osi^Si) operator we are correcting. Using eqs. (|346| , P47|) we get for 
the sum of the four diagrams 
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The diagrams depicted in fig. ^formally describe annihilation matrix elements {H\0{Q)\H) , 
but the result is identical for the production ones {O^ {Q)). Recalling the definition of 0[i^s]{^Pj) 
(see Appendix ^) we can write 



{H\OsCSi)\H) = {H\Os('Si)\H)norn + 
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J=0 J=0 



(148) 



The presence of the ultraviolet singularity indicates that the {H\Os{^Si)\H) operator needs renor- 
malization. The suitable counterterm in the MS scheme is such that the relation between the bare 
(D-dimensional) and the renormalized matrix element reads 
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(149) 



fix being the NRQCD renormalization scale. Since the renormalized matrix elements on the right 
hand side have mass dimension three while the bare D-dimensional one on the left has dimension 
D—1 (a factor of 2{D—1) coming from the four spinor fields of the operator plus a factor of —{D — 1) 
from the nonrelativistic normalization of the \H) states, {H{p)\H{p')) = {2Tr)^~^6^^^{p — p')), 
the /iA~^^ compensates for the difference. 

Using eqs. ( |148| , |149| ) (no renormalization is necessary at this order for the {H\0[i^s]{^Pj)\H) 
matrix elements) we find 
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This equation can now be used to cancel the infrared pole in the ^Pj qqg processes given 
in eqs. (|116| , |120|) : summing to these decays the colour-octet one ^5*^^^ —>■ qq and substituting 
the bare matrix element with the expression given above one finds (restricting ourselves to the 
infrared singular parts): 

j:^^eP?^qq9)+r"CS,^'^^qq) = 



J2 I --UfCpal 



1 MM') fAT^fi'Y r(i-e) 



K 



M2 J r(2-2e) 



+ a J 



+ IR-finite 



itel 



a^TT Mvr/i^Y r(l-e) (1 - e) 



4m2 / r(2 - 2e) (3 - 2e) 



1 



2e 



X + In47r-7B — 

, 9 



4as 



nfCFlJ.^'^{H\OiCPj)\H) + IR-finite 
J 



, 1 , 



+ 



:i5ii 



7 Production processes 

7.1 Kinematics and factorization of soft and collinear singularities 

The technique used to extract the soft and collinear singularities from the real emission processes is 
similar to the one employed in the study of decays. The kinematics of the process ki + k2 P + k^, 
where P is the momentum of the heavy quark pair ( identified by ) and ki are the momenta 
of the massless partons, can be described in terms of the standard Mandelstam variables s, t and 
M, defined by 

s = {ki + k2)\ (152) 
t = {k^-k^f = -■l{l-x){l-y) , (153) 

u = {k,-k,y = -'-{l-x){l + y) . (154) 

Here we also have introduced the Lorentz invariant dimensionless variables x = AirP' / s and y 
(— 1 < ?/ < 1), defined by the above equations. In the center-of-mass frame of the partonic 
collisions, the variable y becomes the cosine of the scattering angle 9. In terms of x and y the 
total partonic cross section can be written in D dimensions as follows: 

^ = J d'^mi^^y) Mx,y) , (155) 

where M. = is the spin- and colour-averaged matrix element squared in D dimensions and 

(i$(2)(a;, y) is the D-dimensional two-body phase space: 

rf$(2)(x,i/) = - (-] T{l + e) — {l-xy-'^{l~y')-^dy. (156) 
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The soft and collinear singularities are associated to the vanishing of t or u, which appear at most 
as single poles in the expression of A4. One can therefore introduce the finite, rescaled amplitude 
squared ^A■. 

ut s'^il — xY{l — y"^) ^ ^ 

In terms of M.^ the partonic cross section reads as follows: 

cr(x) = — / dy{l-y^r^-^M{x,y)dy, (158) 

s-^ J-l 
4' /47r\' 1 

^-If (t) r<^ + ^'32^- <™' 

As in the case of decays, soft and collinear singularities are now all exhibited by the universal 
poles which develop as x ^ 1 and y"^ ^ 1. The residues of these poles can be derived without 
an explicit calculation of the matrix elements, as they only depend on the universal structure of 
collinear and soft singularities. We will carry out an explicit evaluation of these residues in the 
case of colour-singlet production in the gg Q^^^g process in the next subsection. The other 
cases are similar, and will be discussed with fewer details in the following. 

7.1.1 gg gQ^^^ processes 

We start by considering the soft limit, x — 1. The following distributional identity holds for small 
(1 - X)--- = -^6(1 -x)+ (^) - 2e M'-^A + Oie^) (160) 



2e \1 — xj o \1 — X 

where p = M^/Shad, P = " Pi ^^^1 the p-distributions are defined by: 

1 1-1 

dx [d{x)] t{x) = / dxd{x) [t{x) - t{l)] . (161) 
p J p 

We can therefore write, with obvious notation, 

cr(x) = (r^=i + (Ta;^i . (162) 
The first terms on the right-hand side is given by the following expression: 

a.=i = J_^dy{l~yY'-^Mix = l,y), (163) 

The X — 1 limit of A4 can be easily derived from eq.(^): 

M{x,y) ^ g^CA-MBorn (164) 

ut 

M{x,y) ^ 4sg^CAMBorn, (165) 

where A^Bom is the £)-dimensional Born amplitude squared for the gg ^ Q process, which is 
independent of y. The integration over y of eq.( |163| ) is elementary, and leads to the following 
result: 

= IMhCa -P-'' CTBorn , (166) 

35 




ouuuuu 

a 



sJUUUUU 



sJLfiJLSUU 



OLAJLJUU 



sJLaJLSUU 

c 








Figure 6: Diagrams for the real corrections to the gg channels. Per- 
mutations of outgoing gluons and/or reversal of fermion lines are always 
implied. 



where H is defined by 



H 



Til-e] 



TT 



r(i + e)r(i-2e) i-T^^ + f(^'). 



and (TBorn is the D-dimensional Born cross section: 
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The collinear singularities remaining in ax^i can be factored out by using the following distri- 
butional identity: 
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where the distributions on the right-hand side are defined by: 



(169) 



;i7o) 



36 



The contribution Ox^^i can then be spht into three terms: 
The term cTfinite has no residual divergences, and is given by the following expression: 
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(172) 



This piece explicitly depends on the nature of the quarkonium state produced. For processes 
whose Born contribution vanishes, this is the only non zero term. 

The remaining pieces, o"y=±i, are given by: 
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The limits for ?/ — > ±1 of Jv[[x^ y) are universal, thanks to the factorization of collinear singulari- 
ties: 
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where Pgg{x) is the standard D = A Altarelli-Parisi splitting kernel, defined in Appendix Using 
these relations we get: 
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where 



- = - - 7e + log(47r) 
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The collinear poles take the form dictated by the factorization theorem. According to this the 
parton cross section can be written as: 



d0-ij(pi,P2) = J2^^kl{XlPl,X2P2)'i^kiiXi)Tij{x2)dXidX2 , 
k,l 

Tij{x) = 5ij5{l-x) - - ^(-^l Vij{x) + Kij{x) 



(177) 
(178) 



where da is free of collinear singularities as e ^ 0. Here we allowed the factorization scale /ip 
to differ from the renormalization scale fi. The functions Vij{x) are the D = 4 Altarelli-Parisi 
splitting kernels, collected in Appendix ^ and the factors Kij are arbitrary functions, defining 
the factorization scheme. In this paper we adopt the MS factorization, in which Kij{x) = for all 
i,j. For the definition oi Kij{x) in the DIS scheme, see for example ref. 1^^. The collinear factors 
r(a;) are usually reabsorbed into the hadronic parton densities, and the physical cross section is 
then expressed as: 

dcrHrH2iPl,P2) = J2^^kliXiPi,X2P2,^F)FkHiiXi,fXF)FiH^{x2,(J,F)dXidX2, (179) 
k,l 
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where /ip) is the density of the parton k in the hadron H, evaluated at the factorization 

scale /ip. Expanding eq. (|177|) order by order in a^, we extract the following counter-terms cr^=±i: 
defined by 



a. 



y=±l 



1 as I 
e 2tx 



'Pgg{x)xaQ 
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+ (6o + 4CAlog/5)5(l 
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Putting together all pieces, we come to the final result for the real emission cross section: 
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XVgg{x) log ^ + 2X(1 - X)Pgg(X) 



'log(l — x) 
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where Colgg Q'^'] is the D- dimensional, Born-level partonic cross section for the production 
of the quarkonium state H via the Q^^' intermediate state, after removal of the 5(1 — x) term 
(see eq. ( |168|) ). The finite functions fgg{x), obtained from the explicit evaluation of eq. (|172|) , are 
collected in Appendix |C[ 

The result for the ^5*0^^ state agrees with previous calculations |^. The results for the P 
waves are new. In the case of ^Si and ^Pi production the O(a^) process is the leading order one. 
Integration over the emitted gluon phase space is completely finite, and the results have been 
known for long time in the literature (see e.g. ref. 0). For completeness, we collect them in 
Appendix |C[ 



7.1.2 qg qQ^^^ processes 

The Born level processes qg Q'^' identically vanish. As a result no IR divergence is present 
at 0{a^) and virtual corrections are not present. The only singularities appearing at this order 
come from the emission of the final-state quark collinear to the initial-state one. The behaviour 
of the amplitude in the y ^ 1 collinear limit is again controlled by the Altarelli-Parisi splitting 
functions: 

AT(x,y) ^ 8nsa^^^Pgg{x)^^MBorn . (183) 
In analogy to the gg case, one introduces the following counter-term in the MS scheme: 

41 = - ^(^] Vgg{x)xao, (184) 
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where Vgq{x is) defined in Appendix ^ 

Following a procedure analogous to the one detailed in the case of gg production, we find the 
following result: 



X 



|p,,(x)log^^i-2^ 



c, 



X 



[Q 



[1] 



IqW 3p[l] 3p[l]l 
^0 5 -'0 5 -'2 J 



(185) 



where the functions fgq{x) are collected in Appendix y, together with the result for ^Pi^^ pro- 



duction, for which no coUinear singularity is present to start with. Our result for ^5*0^' production 
agrees with those presented in refs. |^. In the case of ^Pq^^ and ^P2^^ we differ from ref. [0 
by the factor Cpx"^ /2 in eg. ( |185| ). This piece arises from the MS factorization prescription, which 
dictates use of the D = A Altarelli-Parisi kernel in the coUinear counter-term, eq. ( |184| ). 



7.1.3 qq 



processes 



In this case no coUinear divergences are present upon integration over the final-parton phase space. 
No virtual corrections are present either. Infrared divergences however appear, associated to the 
presence of an intermediate ^S'f ^ state: 



3p[l] 



g] 



■^'^^^p-''us)a^'^8il-x)ion^p,)) 

Dp 9(2m) Ve 6 J 



(2m) 



1 



1 — X 



WPP]i^){0^i'Pj)) [J = 0,1,2], (186) 



where Dp is defined in Appendix ^ The infrared divergences is absorbed by the inclusion of the 
process, as discussed in Section p. The resulting finite expressions and the functions 



qq 



fqq{x), as well as the trivial result for the ^Sq'^ production, are collected in Appendix 0. 



7.1.4 



gg 



processes 



The diagrams of fig. y have been evaluated in 4-dimensions for colour-octet production in ref. . 
Since the Born cross sections gg — > Q^^^ do not all vanish we expect both infrared and coUinear 
singularities to appear. After subtraction of the coUinear poles in the MS scheme the partonic 
cross sections read 



{2mf 



(187) 



a^[gg^-'PPg] 



(2m)' 



(188) 
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a be 



Figure 7: Diagrams for the gq channels. Reversal of fermion lines is 
always implied. 
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+ 
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XVgg{x) log + 2X{1 - X)Pgg{x) ^ ^^J^^ 



1 — X 
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with the functions fgg[Q^^^]ix) collected in Appendix y. 



;i89) 



7.1.5 qg gQ^®' processes 

These processes proceed via the diagrams of fig. |^. Diagrams b) and c) only contribute to ^Sf^ 
production. 

Since no 0{al) cross sections exist for these channels the results are expected to be infrared 
finite, only collinear singularities being allowed as long as gg Q'^' is non vanishing. Notice the 
special case of ^S'f production, whose collinear singularity is rather subtracted by the qq ^Sf^ 
Born term. 



The cross sections read: 
a^[gq^^S?^q] - 



a 



^^H\„^ _^ 3C'[8] 



-0", 



TT 



sil — x) 



-Vqg{x) log — + Tpx'^il - x) 



+ U?S?]{x)\ (190) 
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Figure 8: Diagrams for the real corrections to the qq channels. Per- 
mutations of outgoing gluons and/or reversal of fermion lines are always 
implied. 
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7.1.6 qq gQ^^^ processes 

The diagrams for these channels are depicted in fig. ^. After subtraction of the coUinear singu- 
larities the cross sections read 



a 



H 



[qq ^ ^Sf^ g] 
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(194) 



and Bp and -D^ are defined in Appendix ^ 
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A few comments are in order. 

Tfie production of ^Sf'^ sliows singularities of both collinear and infrared type. The latter will 
be cancelled by the addition of the virtual corrections to qq ^SfK 

The a^[qq —>■ ^Pf^ g] cross sections only show a singularity which is neither collinear nor 
infrared in the usual sense, i.e. to be cancelled by the addition of virtual corrections. The Born 
cross sections qq — > "^Pf^ are indeed zero, therefore not allowing for such cancellations to take 
place. This singularity is originated by the gluon emission from the heavy quark line, and it can 
only be eliminated, in the spirit of the factorization approach, by its re-absorption into the matrix 
element {Og{^Si)) when adding to this cross section the degenerate one (in the infrared endpoint) 
da"[qq ^^fl] (gee Section 

The process qq —>■ ^So^^^g is completely finite and it is shown in the Appendix ^ together with 
the functions 



7.2 Inclusion of virtual corrections 

In order to obtain the full 0{al) correction we need to calculate the one-loop QCD corrections 



to the four non vanishing Born processes qq — > QQ[^Sf'^], gg QQ[^sll''^\ gg QQI'^Pq''^^] and 
99 ~^ QQl'^P^^'^^]- The relevant diagrams for the qq and gg channels respectively are shown in 
figs. ^ and |^. The diagrams are depicted in fig. H for the qq initiated process and in fig. ^ for the 
gg one. The full virtual correction is constructed from the various parts listed in the tables by: 



TT 
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k 



Explicit expressions for J2k'^kfk can be read out of eqs ( |124| , |128| , 132 ). The sums of real and 
virtual corrections take simple forms. In the case of gluon-initiated colour-singlet production one 
obtains: 



a^igg^Q^'^X] 



cr^[99^Q][S{l-x) + ^ 



B 



3 



LL 2T72 

8Ca log' /? + 2 6o log — + 8Ca log /? log — 



+ 



+ 



4:Tn 
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(196) 



where the coefficients -Bq[i] are the finite parts of the virtual corrections, defined by eq. ( |124| ). A 
similar expression holds for the g g Q'^^ X and qq '^S'f ' X NLO cross sections. The final 
results for the finite sums of real plus virtual corrections are collected in Appendix |C|. 
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8 Conclusions 



We have performed a calculation of next-to-leading-order QCD corrections to total hadronic cross 
sections and to light-liadron decay rates of heavy quarkonium states. Both colour singlet and 
colour octet contributions are included. 

We have extended the technique of covariant projections to D dimensions, to use it within 
dimensional regularization. Where common results exist, they have been found in agreement with 
the ones given by the threshold-expansion technique recently introduced by Braaten and Chen. 

All the singularities which develop during the calculation are found to cancel either via the 
usual QCD mechanisms (soft/virtual cancellation, coUinear factorization, mass and coupling renor- 

malization) or via NRQCD cancellation mechanisms, by taking into account radiative corrections 
to the NRQCD matrix elements. The final formulas are therefore self-contained and do not need 
ad-hoc prescriptions to produce finite results. 

Results for quarkonia photoprodTiction and decay into one photon plus light hadrons can 
be easily extracted from these calculations: explicit results will be presented in a forthcoming 
publication. 

A phenomenological study of quarkonia NLO total production cross sections and decay rates, 
including all the processes calculated in the present work, will be the subject of a separate publi- 
cation. 
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A Symbols and notations 

This Appendix collects the meaning of various symbols which are used throughout the paper. 
Kinematical factors: 



where s is the partonic center of mass energy squared and Shad is the hadronic one. v is the velocity 
of the bound (anti)quark in the quarkonium rest frame, 2v being then the relative velocity of the 
quark and the antiquark. We also define: 



M = 2m , 




(197) 




r(l + e) = 1 + e -7e + log(47r) + log 




(198) 
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and we denote a perturbative QQ state with generic spin and angular momentum quantum num- 
bers and in a colour-singlet or colour-octet state by the symbol 



(199) 



Altarelli-Parisi splitting functions. Several functions related to the AP splitting kernels enter in 
our calculations. We collect here our definitions: 



where 



'Pqqi^ 



Pqg{x 
'Pqgix 



Pgq (-^ 



'Pgqi^ 



Pggi'^ 



Pggi-^ 



Cf 



l-x 
(1 - x), 



-e{l-x) 



Pqq + CF[- + ^\0g(5]5{l-x) 



Tf 
Tf 

Cf 



x^ + xy - 2e x(l - x] 
x^ + (1 - xf^ 



Cf 
2Ca 

:2Ca 



l + {l-xf 

X 

1 + (1 -x)^ 



e X 



X 



X 



.1 — X 
X 



+ h x{l - x) 



X 



{l-x), 



+ 



X 



+ xil — x) 



X 



Pgg{x) + (6o + 4Ca log/3) 5{l - x) 



11 2 
bo = —Ca - -TFTif 



(200) 

(201) 

(202) 

(203) 
(204) 

(205) 

(206) 

(207) 

(208) 
(209) 

(210) 



with Uf number of flavours lighter than the bound one. The Pij are the D-dimensional splitting 
functions which appear in the factorization of coUinear singularities from real emission, while the 
functions Vij are the four- dimensional AP kernels, which enter in the MS coUinear counter-terms. 
The p-distributions are defined by: 



dx [d{x)] t{x) 



dx d{x) [t{x) - t{l)] . 



(211) 



Colour Algebra 



{T",T^} 
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F 



2 



iV2_i _ 4 



Ebc f'^^f^^^ = CaS"^ Ca = N, = 3 (214) 

Df = Y.^n = Nc = 3 (215) 

i 

Da = ^ = A/',2 - 1 = 8 (216) 

a 

The following formulas were found to be useful: 

Y^W^i = l{Su6,,-j^J^,6M) (217) 

Tr(T"r''r^) = ^{d^^'' + ir^') (218) 

Tr{T''{T'',r}) = ^d^^' (219) 

CfDf = Da/2 (220) 

C2{F)Df = BfDa (221) 



NRQCD operators. Although we never make explicit use of them, we collect here for ease of 
reference the definitions of the NRQCD operators relative to the states considered in this work 

ip, 



(see also [15| for details and a comparison with theirs). In the case of 5* states we have: 

1 1 



and for P states: 



Oi('Po 



Oi('5i 
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■■x-x 
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01('P1) = V'^ 



i D. 1 1 t/ ^ D^ 1 1 

=x ■ x'(- 



2 'V2v^ 



3^ ^ 2 V2^ViVe 



XX' 



2 'V2v^ 

i D cr . 1 
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(222) 
(223) 
(224) 
(225) 



(226) 
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= X ■ X'i-- X (228) 

OsCP^) = ^^{--f)^V2T'^x ■ xK--f)^V2T'^i^ (230) 
Osi'Po) = ^^n-^ ■ ^)V2T^XXH-^ ■ (231) 
C^sm) = ^^n-^ X ^)V2T^X ■ XK-^ X ^) V2T"^ (232) 

C?8(^P2) = ^t(-i±L^) V2T"^^t(_i|L^) V2T> (233) 

Writing all previous operators as products of fermion bilinears, 0{n) = 02{n)02{n), we also recall 
the definition of the various NRQCD matrix-elements appearing in the production cross sections: 

(^mH) = j:{0\O^^^^\HX){HX\Oi''>\0). (234) 

X 

Notice that, according to the discussion in Section our conventions differ slightly from those 
introduced in ref. (and labelled here as BBL): 

^1 = i^^r"' (235) 



2 AT, 
2 



Cg = -2C^^L = Of^ . (236) 



B Summary of 0(aj) Results 

The D-dimensional cross sections and the decay rates read 

a{^J - ^^-^l/^'' ^H) = a{^J - ^--^L;^'^')^^^!!^, (237) 

T{H ^ ,_^-) ^ f(2^+iL;^''l ^ zj){H\0[^,s]e'^'Lj)\H) , (238) 

the short distance coefficients a and F having been calculated according to the rules of Section ^ 
Ncoi and Npoi refer to the number of colours and polarization states of the QQ[^^'^^L^j'^^] pair 
produced. They are given by 1 for singlet states or Da = — 1 for octet states, and by the D- 
dimensional A^j's defined in Section ^. Recall that the matrix elements appearing in the equations 
above are meant to be the bare D-dimensional ones. Making use of their correct mass-dimension, 
3 — 2e and 5 — 2e for S and P wave states respectively (see Section]^), gives the right dimensionality 
to D-dimensional cross sections and widths, i.e. 2 — D = —2 + 2e and one. 
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B.l 0{a^) decay rates 

We shall use the short-hand notation 

YH^q[i,8] ^ ^ ^ Q[i,8] _^ (239) 

to indicate the decay of the physical quarkonium state H via the intermediate QQ state = 

The D-dimensional Born level decay rates read: 

r^ornCSo^'^ - 99) = (1 _ _ 2e){H\0^{}S,)\H) (240) 

m 

rforn('^i™ ^99) = ^ (241) 

rgo^n^i-o'^' ^ ^^) = ^^^'"l^'''' ^m^{H\0,{'Po)\H) (242) 

7/7' O 

T^^^A'Pi'^ ^99)=^ (243) 

rL.(^P.^^' ^ ..) - ^^^(^) (L'^lxs^'-t) (244) 

rforn('5i^'^ - = nj^-^^^,^^^{H\0,{'S{)\H) (245) 

rforn('^o^'^ - ^7^7) = 5^^^^J^$(2)(1 - e)(l - 2e){HWS,)\H) (246) 

m 

r^o.n(^i^o''' - 99) = B, ^^^"^f $(.) ^ 08(^Po) (247) 

rforn('A'''-^7^?)= (248) 

where $(2) is the two-body phase space given in the equation (p7|). 



B.2 0{al) cross sections 

We shall use the short-hand notation 

a^(zj ^ Q[i'«]) = (t(«j ^ if) (250) 

to indicate the production process of the physical quarkonium state H via the intermediate QQ 

state Q}^''-^=QW'^^l}^\ 

The D-dimensional Born cross sections read: 

^L..{99 - 'S^') = ^^^^^5(1 - x){0^{^S,)) (251) 

777- 1 — c 
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C Summary of 0(ag) Results 
C.l Decay 



(252) 
(253) 

(254) 
(255) 

(256) 
(257) 

(258) 

(259) 
(260) 



We distinguish the two different cases in which the 0(q;^) subprocess is cither Q ^ qq or Q ^ gg. 
In the former case the inclusive 0(0;^) inclusive annihilation rate into light hadrons (LH) of the 
quarkonium state H through the component Q can be written as follows: 

r^(Q ^ LH) = rf„r,(Q ^ gg) + r"{Q ^ gg"^) + r^(Q ^ ggg) + r^(Q ^ gg) (261) 

In the latter case one has instead: 

r^(Q ^ LH) = r|,,„(Q ^ gg) + T^{Q ^ qqg) + T^{Q ^ ggg) + T^{Q ^ gg) (262) 

Given a quarkonium state H, we remark that, in order to obtain an infrared finite result, it is 
necessary to sum over all the intermediate configurations Q which contribute at any given order 
in the velocity v. Typical cases are given by the '^P/^'- '^5*^' and ^Pf \ ^Sf^ configuration pairs that 
appear at the same order in v for any quarkonium state H. The corresponding NRQCD operators 
then mix under rcnormalization group equation. 

The inclusive decay rate is thus given by 

r{H ^ LH) = ^ r^(Q ^ LH) (263) 
Q 

and we list here the most interesting 0{al) contributions to the annihilation of a state H into 
light hadrons. 
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r^(^^o^ ^ ^ LH) 

r^(i^o[«] ^ LH) 

r^(3^jii ^ LH) 

r^(35jsi ^ LH) 
r^(3p^W ^ LH) 

r^(3p^W ^ LH) 
r^(3p2W ^ LH) 
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We remind the reader that every time the NRQCD factorization scale fi\ appears the r^('^S'i'^' 
LH) channel has to be added for consistency, and the {H\Os{^Si)\H)^'^'^^ matrix element is then 
understood to be the renormahzed one (see Section ^. 

We make now some comments on the comparison of our results with previous calculations. The 
calculations for ^Sg^' decays have already been performed in refs. [|, Our results agree with 



theirs. The calculation of the 0{a^) decay rate for ^Pq and first appeared in ref. 0, where 
a massive-gluon regulator was used to treat the IR divergencies. Ours is the first independent 
calculation since then. The results differ in the Ca terms, for which we find: 



where Barbieri et al. find instead: 
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The differences amount to a factor Ca/3 and Ca/4 for J = and J = 2, respectively. We point 
out tliat tlie 0{as) corrections to tlie two-plioton decay rates of '^Pq^' and ^^2^', wfiich can be 
obtained by setting Ca = and by summing diagrams a tlirough d of the virtual corrections 
displayed in Tables |^ and |], coincide with the results of ref . . 

The calculation of the 0{a^) decay rate for ^5*0^^ has already been performed in ref. |^0[. The 
results differ in the terms, for which we find: 

while Huang et al. find: 

which differs by a factor of C^/A from ours. A comparison with the diagram-by-diagram break- 
down of the calculation performed in ref. shows agreement in the evaluation of the virtual 



corrections, and the discrepancy therefore originates from the real-emission partQ. 

C.2 Production 

We define: 

cr^i^J - Qm -x)^ ai,,S3 Q) (280) 
The 0{a^g) cross sections are given as a function of the variable x = M^/s. 

C.2.1 Colour-singlet channels 
The gg — > Q^X channels 

-T^bs^'sWal = ^^„?S?]^x){0»eS^)) (281) 
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D\ 3(-l +x)2(l +x)3 
+2x^(5 + 2x + x2) log xl (282) 



^""[gg^'P^g] = ^WP?]{^){on'Pi)) (283) 



„ rsoWu ^ CaCf 256 
fa9[Pi\{x) = 



10 - 8x - 157x2 - 9x^ - + 136x^ 



D\ 9(-l +x)4(l + x)5 
-F404x^ + 26x^ - 212x^ - 144x^ - 39x^° - x^^ + (-48^^ 
-48x^ - 480x^ + 84x^ + 204x'^ + 336x^ + 312x^ 

+ 108x^ + 12x^°)logx] (284) 



^After this paper was released as a preprint, the authors of ref. |Q reviewed their calculation. We have been 
informed that their final result now coincides with ours. 
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Where: 



AtotpPi''] 



+ 



+ 



log ^ + 2a;(l - x)P,,{x) f 



1 -X 



[QW=i4'^,=^pW,3pW] (285) 



-5 



2m 



+2 6o log — + 8Ca log /3 log — + SC^ log" (3 

/ip /ip 



2m 



+2 6o log ^ + SCa log /? log — + SCa log' /? 

//p 

-4C. + C.(i-^ + ^log2) 
+2 6o log + 8Ca log /3 log — + 8C^ log' /3 



6(-l + 

+9x^ - llx^ + 12(-1 + 5x' + 2^3 
+x^ + 3x^ + 2a;^)logx) 



12 + llx' + 24x^ - 21x^ - 2Ax^ 



(286) 



(287) 



(288) 



(289) 



172 - 56x - 789x' + 244x^ + 1721x^ 

54(-l + x)^(l +x)5 
-696x^ - 802x^ + 560a;^ - 210x^ + 80x'^ + 7x^° - 1322;^^ 
-99a;^' + 12(-9 + 31x' - 14x^ - 40x^ + lOx^ + 176x^ 
-42x^ - Ix^ - IQx^ + 41x^° + 24x^^) log a; ] 



(290) 



c. 



106 - 32x - 207x' + 271x^ + 752a;^ 



36(-l + a;)3(l +a;)5 
-1032a;^ - 256a;^ + 266a;^ - 360a;^ + 728a;^ + 31x^° - 201x^^ 
-66x^' + 12(-6 + 22a;' - 8x^ + 74a;^ + 31x^ + llx^ 



-204x^ + 86a;^ + 17a;^ + 5a;^° + 12a;^^) log a; 
The gq{q) Q^^^X channels 



(291) 



a. 



TT 
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where 



X 



Pggi^) log 



4nr{l — x)' 



X 



X 



(2m)' 



12 



-(-1 + a;)(43 - 14a; + 4x^) + (9 - 9x + Ax'^)\ogx 



[-1 + x) (53 - 16x + 20x^) + (12 - 12x + 5x^) log a; 



Cf 128 
DfDa 3 



1^-1 + a;) (-5 + 4a; + 4a;^) - ar^loga; 



The qq — > Q^^lX channels 



(292) 
(293) 

(294) 

(295) 
(296) 

(297) 

(298) 



where 



'P?9] 



(2m)5 



+ 



g-(,,.,,_^_,og|^).(l-.) 



1 



.1 — a; 
[J = 0,1, 2] 



{on'pj)) 



Cf 32 2 



fJ'P 





X 




3 




128 




9 


Cf 


256 


Dl 


9 


Cf 


256 



Dl 45 



a;^(l -x) 
x^{l - 2>xf 
x^{l + x) 
a;2(l + 3a; + 6a;^) 



(299) 
(300) 



(301) 

(302) 
(303) 
(304) 
(305) 
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Note that, alike the T^ppj^' qqg] process, a^[qq ^Pj g] is only infrared finite when at 
least the Born term of the cr^[qq — > ^-Sf] cross section is added to it. We have performed the 
cancellation and an explicit dependence on the NRQCD rcnormalization //a results. The NRQCD 
matrix element in a^[qq ^Si'^] is then understood to be the renormalized, //a scale dependent 
one 



C.2.2 Colour-octet channels 
The gg Q)-^^X channels 



a^igg^^Pl'^g] 



(2m)" 



(306) 



— — \:r- [l08 + 153x + AOOx^ + 65x^ 

36(-l + x)2(l + x)3 L 

- 356x^ - lS9x^ - 152x^ - 29x^ + (l08x + 756x^ 



+ 



432x3 + 704x^ + 260x^ + 76x^) log 



X 



(2m) 
CaBf 



(307) 



(30J 



256 



Da 



9 } 777 7T 10 - 8x - 165x^ - - 182x^ + 220x^ 

m2x^ + 44a;^ - 276a;^ - 244a;^ - 79x^° - Ax^^ + (-48a;2 - ASx^ 

:^ + 144a;^ + 48a;'^ + 576a;^ + 540a;^ + 192a;^ + 24x^°) logx] (309) 



- 660x^ + 



a^igg^Q^'^X] = a^[gg ^ Qy^\]i^5{l - x) + {KUQ'^'^K^ - x) 



+ 



1 -X 



(310) 



where: 



A [3 pi' 
^tot[ -r 



3p[8]l 



+2 6o log — + 8Ca log /5 log 2Ca log P + 8Ca log' /? 

^ / 7 7r2\ ^ /71 35, \ 

^M-3 + T +^M 54 - 24 + 27^°S2 



(311) 



54 



■3p[8]l 



2m 



+2 bo log ^ + 8Ca log plog 2Ca log /3 + SCa log" /3 

+ - + - + -log2 



2m 



+2 6o log — + SCa log /3 log 2Ca log /3 + 8Ca log" /3 



12 + 23x^ + 30x^ -21x^- 2Ax^ + 9x^ - 



(312) 



(313) 



+ x){l + 

- 23x^ + (-12 + 6Qx^ + 24x^ + 36x^ + 60a;'' + 24x'^) logx 

— ^4t7 ^ [172 - 56x - 785x2 ^ 254x3 + ISSlx^ - 948x^ 

5A{-1 + x)^{l + xf L 

- 334x^ + 632a;^ - 790x^ + 268x^ + 63x^° - 150x^^ - 207^^^ + (-108 
+ 372^2 - 168x3 - 540x^ - 180x^ + 3348x^ - 804x^ + 312x^ - 252x^ 



(314) 



+ 840x^° + 252x^^) logx' 



(315) 



36(-l + x)3(l +x)5 
- 400x^ - 76x^ - 424x^ + 1252x^ - 153x 



106 - 32x - 215x2 + 260x3 + 1224x^ - 1284x^ 



10 



120,T 



11 



138x^2 ^ 



+ 264x2 - 96x3 ^ io80x^ + 36x^ + 1404x^ - 3180x^ + 1320x^ - 612x^ 



+ 516x'" + 108x'')logx 



(316) 



The gq{q) Q^^^X channels 



TT 



X 



X 



4m2(l - x)^ 



P,,(x) log I] + T^x2(l - x) 



+ Wsf]{x)\ (317) 



a"[gq^^Pf\] 



3^2 



(2m) 



(318) 



aH[gq^Qn q] 



TT 

[Q[«l = ^5fl,3p[«l 3p[8l] 



X „ , , , 4m2(l - x)2 ^ 
-P,,(x)log \-^ + Cf 



X 



X/Xp 



(319) 
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f r3p[8]i _ ^ 



-^(l + 3x)(-l+x)x 



CaDf 



Da 

Bp m 

DfDa^ 



-1 + x){2 + X + 2x^) — x{l + x) log X 
(-1 + x)(-5 + 4x + 4x^) - log X 



CF(-l + a;)(l-loga;) 

Cf_ n 

9" L3 

Cf n 



12 



1 + x) (43 - lAx + 4x^) + (9 - 9x + 4a;^) log x 
1 + a;)(53 - 16a; + 2Qx^) + (12 - 12a; + bx^) log a; 



Q^^^X channels 



- 'M'^] (^(1 + 



TT 



-10 n/Tp 



9 



+Ca y-- + -log2-2 1og/3 



TT 



2m 



+Cf -8 + — + 8 log^/3 + (3 + 81og/3) log — 



+ 2 6o log 



2m 



5(1 -x) 



+ 



+ 



xP,,{x) log ^ + C^a:(l - a:) + 2a;(l - a;)P,5(a:) f 

XjJ,^ \ 1 



1-x 



a 



(2m)^ 
+ 



^^(21og/3---log-)5(l-a;) 



1 



a; 



[J = 0,1, 2] 



X ricMi -Sf 32 X 
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/,_[35f]] = -Cfx{1 - xf - + x + x^) (329) 

U'P'^'] = ^^^\l-^^? (330) 

/,.N^'] = ||^x^(l + a:) (331) 

/..P^f] = ||^^^(l + 3x + 6x^) (332) 

In a way identical to the colour singlet case, [qq '^Pf^ g\ is only infrared finite when at 
least the Born term of the a^[qq ^Sf\ cross section is added to it. We have performed the 
cancellation and an explicit dependence on the NRQCD renormalization /xa results. The NRQCD 
matrix element in a^[qq "^Sf"^] is then understood to be the renormalized, /iA scale dependent 
one (Of (35i))(^'^). 



D Virtual contributions with the threshold expansion 



In this Appendix , we show how the threshold expansion method ^ is used to calculate the 
imaginary parts in the short-distance coefficients for the virtual contributions at NLO in as- Being 
the S'-wave calculations analogous but much simpler, we will sketch here only the derivation for 
P-waves. We start considering QQ scattering in an arbitrary frame in which the pair has total 
momentum P and the relative momentum of the two quarks is small compared to the heavy quark 
mass m. The momenta p and p of the Q and Q can be written as 



p = iP + Lq , 
p = - Lq 



The components of the momenta P and Lq in the CM frame of the pair are 



CM 



CM 



{2E„ 0) 
(0, q) . 



(333) 
(334) 



(335) 
(336) 



When boosted to an arbitrary frame in which the pair has total spacial momentum P, these 
momenta are 



P^ = [^AE^ + P2, P 



(337) 
(338) 



The boost matrix L'^-, which has one Lorentz index and one Cartesian index, has components 



2E„ 



pi pj PO pi pj 

5'^ \ . 

P2 ^ 2Eq P2 



(339) 
(340) 
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These expressions are of course valid only for boosts with a vanishing time component. We first 
calculate the Born amplitude in D-dimensions for the process QQ gg, i.e. diagrams Di,D2 of 
fig. It reads ( all the momenta are outgoing ) : 



QQ^gg 



gl^?' v{-p) 



brpa 



2p-k2 



+ 



Y 



+ m)7'' 



rj-iarj-ib 



u{-p) 



(341) 



We then proceed by two main steps. We consider the following expressions of the spinors of 
the QQ pair in the CMS , 



'E + m 
2E 



lE + m 
2E 



( 



qcr 



E+m I 



V 



7] 



(342) 



(343) 



where ^ and t] are 2-component spinors with suppressed colour indices. Once boosted to a frame 
where the pair has total 3- momentum P, it is straightforward to obtain the relevant independent 
quantities up to the second order in v , that can be formed by sandwiching up to three Dirac 
matrices between v{—p) and u{—p) (v = q/i? and v' = q'/E). They read 



v{—p)u{—p) 
v{—p)Yu{—p) 

v{-P)iYY-YrM-p) 



v{-p){YYl^ - 7V7'')m(-p) 



A 1^ /XN 



1 

m 

X 

1 

m 

X 



1 - v^) ri^a^i 



2 

P^'L'^X^^ + L^^L^.P^ + L^^P^L' 



(344) 



We note that these formulas are basically the ones given in Appendix A of Ref. except that 
only the terms that are specific to the case of P- waves, are kept. The normalization for the 
spinors is however the non-relativistic one adopted in Ref. [15|. Once the amplitude in (|341|) is 
simplified and rewritten using the above amplitudes, one is ready to perform the second step, i.e. 
to calculate: 



where ki and k2 are the momenta of the gluons, Jf^ 



Ik — I (27r)^2fco 

dimensional integral over the phase space associated with the momentum k, the sum is over their 
spins and colours, and the factor of ^ comes from Bose statistics. 



denotes the = (£> - 1) 
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The colour structure is easily simplified by 

rparpb ^ rpbrpa 



fiV^ - 1) iV^ - 2 

^ ^ 1 ® 1 + T'^ (g) 



rj-iarj-ib ^ rj-iarj-ib 



(N! - 1) 



1® 1 



while the phase space integration is performed using the formulas 



-$f2) 



(D-l 



$(2); 



where 



$(2) = — 



1 f A-K \— Ti 



D-2^ 



Svr Um2/ r(L)-2) 



(346) 
(347) 

(348) 
(349) 

(350) 



is the the phase space of two massless particles in D dimensions. In so doing, we eventually obtain 
the following expression: 



ImA^ 



{Nl - 1) 



(D-1)(D + 1) " 



X 



(5D + 2) ^'V • CTT^'r/V ■ (T^ + {2D' - 3D - 8) V ■ v^'Vr^'r^V^ 



+ (2^2 - 3D - 8)^'W ■ <tV^/^v' • 



(351) 



Next we isolate the individual P-wave contributions. To this aim we follow the same procedure 



of section II and Appendix A of Ref. |T5[, but generalizing it to arbitrary dimensions. This can be 
accomplished by first noting that any direct product of cartesian vectors in spatial dimensions 
may be written as 



{D - If 



a ■ b 

D-l 

1 



5'^ + 



{D-2){D + l) 



a ■ b 

D-l 



+ {a^ -a^b')/2 
{D-l){D-2) 



(352) 
(353) 



Applying the above decomposition to the case at hand, we obtain from the NRQCD lagrangian 
the relevant coefficients for the ^Pj states in D — 1 spatial dimensions 



2Nr, 



D- 1 



2A', 



+ 



fl{'P2)-fli'Pl) V-0-®v'-0- 



2A', 



(354) 
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Comparing this relation to ( p51| ) we can read off the imaginary parts in = 4 — 2e dimensions 



Im/i(3Po)o 
Im/i(=^Pi)o 



9(1 -e) 
3-2e 



2 ,.4e 



-$(2) a,fi 



Im/i(3p, 



2)0 



6 - 13e + 4e2 
(3-2e)(5-2e) 



2,.4e 



$(2) a,fi' 



(355) 
(356) 
(357) 



Needless to say, their octet counterparts differ only by the colour coefficient. Once the leading 
order expressions are available, we are ready to perform the next order corrections to (|341|) . We 
can write 

f3: 



lmf,CPj)=lmh('Pj)o 



IT 



(358) 



where the Vj are the sum of all the virtual contributions shown in fig. |^ 

Let us consider, for example, the contribution of graph (e) in fig. ^ In this case the virtual 
amplitude reads: 



A" 



4 4e 



(27r 



D 



-p) Y- 



h + P +rn Q +m 



-r 



^ rj-ibrj-idrj-iC 



F^fAh,-p-Q-h,p + Q] 
{Q + pYiQ + P + hy 



(359) 



where Q is the loop-momentum and F^^'^ is the usual three gluon vertex. By making use of the 
Dirac equation and contracting all Lorentz indices, it is tedious but straightforward to rewrite 
(|359|) in terms of the basic amplitudes in (|344|) . Due to the great number of terms which appear 
in these expressions, all the algebraic manipulations have been performed using MATHEMATICA 
and the package FEYNCALC Only at this stage, loop tensor integrals of the type 



(27r)^ (g2 _ ml) (g + p)2(g + p + 



(360) 



are decomposed with the usual Passarino-Veltman technique, which is performed by computer. 
At the end of this procedure one is left with a tensor structure based on the external momenta 
and a small number of scalar integrals of bosonic type. We evaluated these integrals using the 
results from ref. ||50|. Once all loop integrals are calculated, the remaining internal Lorentz indices 
are contracted and the phase space integration 



virt 

QQ^gg 



(361) 



performed. The colour structure globally factorizes both in the singlet and the octet contribu- 
tions. The final results coincide diagram-by-diagram with those obtained using the D-dimensional 
covariant-projector technique, and are shown in Tables |] and ||. 
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E Fragmentation functions 



We consider here another apphcation of the results of Section ^ We shall deal with the cancel- 
lation of the infrared singularity in the gluon to xj states fragmentation functions, which have 
already been considered in literature in various occasions and in one case fully worked out 

within the matching method with dimensional regularization |^2|. This is a fourth independent 



calculation: our result coincides with ref. [H2 



The general expression for this fragmentation function needs for consistency both a singlet 
and an octet part. Calculations made with the projection method and dimensional regularization 
according to the rules given in Section ^ give the following result in D = 4 — 2e dimensions: 



27 



1 



+ 



2z 



{OrCPj)) 



26 



8(3 - 2e)m- 



Ml-z){OrCS^)) 



(362) 
(363) 



The Pj{z) are finite functions which are of no concern to us here, and can be found in the 
literature, for instance in [^. Nj is the number of polarization states of ^Pj, in four dimensions. 
The constants aj take the values 

. _ 1 . _ 5 _ 19 

- ~6' - "12' - "60- 

The short distance coefficient of the colour singlet matrix element can be seen to be infrared sin- 
gular. This singularity is cancelled by substituting, according to eq. (|150|) , the bare D-dimensional 
{Og ' {^Si)) matrix element which appears above. The final result then reads 



(364) 



with 



d({z,fiA) 



dfiiz) 



z) = d(iz,^A){OrCPj)) + dsiz){orcs^)y>^-^ 



27m^ 



-In^) 5(1 



2m J 



-z) + 



(l-z). 



was 
24m3 



6(1 



and 



«1 = o> 



"^=40- 



(365) 

(366) 
(367) 

(368) 



Finally, differentiating (|365|) with respect to In/XA returns the evolution equation for the 
(0^-'(^S'i)V''a) matrix element: 



d 



djjj 



32xXs(jj_a) 



{orcpj)). 



(369) 



This corresponds to the standard result |jT5|, once the difference in the normalization of the 
NRQCD operators, eq. ( |235| ), is taken into account. 
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